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Refractive Index of Cesium Bromide for Ultraviolet, 


Visible, and Infrared Wavelengths 


William S. Rodney and Robert J. Spindler 


The index of refraction of cesium bromide was measured at 37 wavelengths from 0.365 


to 39.22 microns 
covered with a single instrument. 
decimal place over the wavelength range. 
KRS-5 beyond 20 microns; 


The minimum deviation method was used, 


and the entire range was 


The index changes approximately two units in the first 
The dispersion compares favorably with that of 
and when the effects of inhomogeneity 


and reflection losses are 


considered, the resolving power of CsBr is probably better. 


1. Introduction 


Larg¢ crystals of cesium bromide of reasonably 
od optical quality have recently been successfully 
“own pam xe a new material for infrared studies 
the range beyond the 25-micron (x) limit of KBr 
d out to about 40 u, wherein lie many of the funda- 
ental modes of vibrations of molecules. A mixed 
vstal of thallium bromide-iodide, known as KRS-5, 
as previously the only material available for use 
this region. 
In order to utilize fully any dispersive medium, 
bectroscopists must have a knowledge of the indices 
refraction and dispersion for all wavelengths 
ansmitted by the medium. Such data are also 
eful to physicists for evaluating theoretical dis- 
rsion equations and for studying the forces between 
e constituents of the erystal. The alkali-halides, 
ving the cubic structure, are favorable subjects 
r such studies. 
The authors are fortunate in having access to two 
mples of cesium bromide whose faces are about 
sq in. One of these samples was grown by the 
arshaw Chemical Co. of Cleveland, Ohio, and the 
her was grown at the National Bureau of Standards 
Francis P. Phelps of the Mineral Products Divi- 
m. The refractive indices of each of these samples 
ere determined for 37 wavelengths ranging from 
the latter being near the infrared 
ansmission limit of cesium bromide. 


365 to 39.22 yu, 


2. Instrument 


The instrument (fig. 1) used in these experiments 
a Gaertner precision spectrometer adapted for 
easuring indices of refraction for nonvisible radia- 
m. The telescope and collimator objectives have 
en replaced by mirrors of the same focal length. 
he telescope eyepiece is replaced by a second, or 
it, slit. The radiation is focussed on the exit slit 
; the telescope mirror, and the image of the slit 
the prism face is formed on the detector by use of 
other mirror or a lens of KRS-5. 

Che infrared detectors employed are the lead sulfide 
otoconducting cell for the near ultraviolet, the 


visible, and the infrared to about 2.5 uw; and the 
Golay pneumatic detector beyond 2.5 uw to the limit 
of transmission of the window employed with this 
cell, which in this case is about 40 u. The signal, 
chopped at the rate of 10 times a second, is amplified 
by a gated amplifier controlled by a photocell in 
the chopper unit. The amplified signal is recorded 
on a recording potentiometer. 

The spectrometer is equipped with a set of gears, 
also shown in figure 1, in such ratio that when engaged 
the prism table rotates at one-half the rotation rate 
of the telescope and microscope ring. The gears 
may be readily engaged at will, and they maintain a 
condition of minimum deviation once it has been 
established for any line. The method of minimum 
deviation provides desirable features of high accuracy 
and simplicity of calculation as compared, for example, 
with methods where a constant angle of incidence is 
employed. 


Ficure 1. 


Gaertner precision spectrometer 

4, Auxiliary telescope is used to level and center the prism; B, gear system is 
used to maintain minimum deviation; C, driving mechanism for scanning 
Entrance slit is partially hidden by telescope directly above driving 


mechanism. Collimator mirror is hidden by the telescope. Exit slit 
is seen between prism and auxilary telescope 


mirror 


vuxiliary 
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3. Procedure 


To begin a series of measurements the refracting 
angle of the prism and the deviation angles at 
minimum deviation for several visible lines are 
measured on a Watts precision spectrometer by the 
usual methods and described in a previous paper [1].' 
Index of refraction values accurate to +1X10~° are 
easily attainable by this method [2]. Then an 
auxiliary telescope, figure 1, is used to level and 
center the prism with respect to the entle ‘al axis of 
the mirror system on the Gaertner spectrometer. 
Another auxiliary telescope is used to set a given 
line, usually the 0.6438 u line of cadmium, at its 
minimum deviation position. The table is now 
clamped in this position. The mirror, acting as the 
telescope objective, is brought into the beam so that 
an image of this line falls on the exit slit, causing a 
deflection of the potentiometer pen. The scale 
position corresponding to maximum deflection of 
the pen is observed with the microscopes. This is 
repeated several times, and the microscope is set at 
the average of these readings and clamped. The 
gears are now engaged and the telescope and prism- 
table clamps released. A condition of minimum 
deviation will then prevail for other spectral lines. 

Unfortunately, with this mirror system, one 
cannot directly observe the position of an undevi- 
ated beam or measure twice the minimum deviation, 
as is feasible in visual measurements. It is therefore 
necessary to compute the average reading for an 
undeviated beam by applying a few deviations, as 
observed visually on the Watts instrument, to the 
scale readings on the Gaertner spectrometer for the 
corresponding lines of the visual spectrum as deter- 
mined by the use of the detector. 

The spectra are now scanned by using a driving 
mechanism consisting of a synchronous motor and a 
gear segment attached to the telescope assembly. 
The spectra used are the emission lines of mercury 
and cadmium for the ultraviolet, visible, and 
infrared to about 2.3 yu, the absorption bands of 
polystyrene from approximately 3 to 15 yu [3], 1,2,4- 
trichlorobenzene from approximately 15 to 20 u [3], 
and water-vapor bands for the rest of the range [4]. 
The bands of carbon dioxide at 4.2 and 14.9 yu are 
also used. The scanning serves to identify and 


locate the approximate scale position corresponding | 


to these lines and bands. The actual measurements 
are made by the method used to determine precisely 
the seale position for the 0.6438 u line. For some 
broad bands and for the region of low inte nsity, scale 
positions may be read from the graphs. This is 
done by using a relay to mark the graph at intervals 


varying from 1 min of are to 15 see of arc, depending | 


on the speed at which the telescope is being driven. 


4. Data 


The refractive indices of the two crystals of cesium 
bromide were measured at room temperatures near 
24° and 31° C. The temperature was determined 
———$—————— 

! Figures in brackets indicate the literature references at the end of this paper. 
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by placing a calibrated mercury 


ture at 1-min intervals with a telescope and ayo, 
ing. This procedure was repeated sever; 


each day, and the variations between these ima 
Indices were det 


did not exceed +0.2 deg C. 
mined for 37 wavelengths, ranging from 0.36; 
39.22 yw, at each temperature. Temperature 


efficients of refractive index were determined. y 
the indices at 24° and 31° C were adjusted to x 
and averaged for each sample. These ayen 


values were again averaged, giving the averag; 
both samples. 


for both samples is 7.9 10~*/deg C. 

Table 1 lists the indic es at “this temperatur 
obtained from the observations and as comput 
by means of eq 1. 


:, :' Ts 2 
2 2 ky? ae 
it Go = Tat 


M m 
ee eo 


where 


$1110.49 
119.96)? 

0.0290764 

(0.15800)2 


\= wavelength, microns M 
a? = 5.640752 L 
k = 0.000003338 

p=0.0018612 


14390.4 


0.024964 


Although this equation has 7 constants, only 


them were determined by means of a simultaneoy 
The constants appearing in the deno: 


solution. 


Observed and com puted data on index of 
of CsBr (27° C) 


TABLE 1 


Observed Computed 
index, index, 
n 


Wavelength 
microns 


OSOLS 75118 
S288 75050 
HH46.56 73344 
435835 72333 
‘M40074 T7OLSY 70193 


643847 69202 69204 
01398 67766 67706 
12st) 75 67584 
52952 7237 67238 


7011 715 67157 


w10 fOs0H) 6886) 
258 HHTU4 H67W3 
465 66587 66620 

9. 724 (6283 66285 
66118 66115 


65504 65599 
65474 65472 
65375 65375 
64967 64974 
64795 64801 


64184 64197 
63846 63856 
63565 63573 
632354 63231 
62817 62802 


. 62521 62520 
62284 62282 
61034 61028 
60749 60713 
60591 60550 


60198 60192 
59584 59576 
S5R835 5SS831 
. 4 58317 
58069 58071 


. 57183 57147 
55990 56106 








the -TTLOMete 
directly over the —_ and observing the | mpen. 


The average temperature coefficiny 
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of two terms represent the infrared and 
et absorption bands. The ultraviolet term, 
letermined by taking a weighted mean of 
measured bands. The infrared term, LZ’, 
mate based on the measured and computed 
ff the infrared absorption bands of cesium 
and is probably low. 
obvious anomalies occur in the residuals, for 
at wavelengths 6.465 and 35.45 uw. It is 
hat the first of these bands consists of several 
ed bands, indicating that the anomalous resid- 
be due to an inaccurate value for the wave- 
\ value of 6.695 uw arrived at by taking a 
tly weighted average gives a more consistent 
The bands at 35.45 and 39.22 yu are quite 
ausing considerable difficulty in locating their 


of refraction as a function of wavelength is 
in figure 2. This graph is particularly inter- 
is it shows the smoothness that results from 
data over a large range on a single instrument 
ngle sample and under the same conditions. 


is hoped that consistent data such as these will 


mselves more readily to theoretical analysis. 


ital change in index amounts to approximately 
nits in the first decimal place over the wavelength 


onsidered. The values plotted are the indices 
Differences between samples are small. 
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Index of refraction of CsBr as a function of wave- 


length. 


Wavelength scale is logarithmically graduated. 


The index values of the Phelps crystal are lower by 

| several units in the fifth decimal than corresponding 
values for the Harshaw crystal. The only prior data 
available are 3 values in the visible spectral region 
given by Sprockoff [5] to only 4 decimal places at an 
unspecified temperature. 

The preliminary values published by the authors 
[6] for only one of the crystals are in substantial agree- 
ment with values published in this paper. There 
are some changes and additions beyond 30 yu, and in 
many cases another significant figure has been added. 

Equation (1) has been used for computing values 
of n—1 for cesium bromide that are listed in table 2. 
They are considered as the best values obtainable 
from the measurements here described and are prob- 
ably accurate to within +1 or 2X10 except for 
wavelengths longer than about 30 x. 

The dispersion An/Ad of CsBr is shown in figure 3. 
Both abscissa and ordinate are logarithmically 
graduated. The dispersion in the far infrared 
increases but remains less than in the visible region 
by a factor of 10. The values of An/AX are plotted 
at the midpoint of the AX increment. 

Values for the dispersion of KBr and KRS—5 
[7, 8] were also computed at various wavelengths 
and plotted for comparison purposes. We see that 
KRS-5 has a higher dispersion than CsBr throughout 
the spectral range considered. KBr has a higher 
dispersion than either from 10 to 25 yu, where it be- 
comes opaque. This result would lead to the con- 
clusion that KRS-5 is better than CsBr as a dis- 
persive material, if the effects of its optical inhomo- 
geneity and larger index could be ignored. 
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Ficure 3. Dispersion as a function ef wavelength. 





Both abscissa and ordinate are logarithmically graduated. 
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Tasie 2. Refractivity, (n-1) 


This table gives the average refrac 


Refractivity 
Wave- 
length 


The results of practical tests on these materials [9] 
indicate that CeBr prisms gives better resolving 
This better resolving power could, largely, 
be due to better optical homogeneity of CsBr. The 
effects of higher index of KRS-—5 also deserve some 
consideration because there are practical limits to 
the size of the crystals grown. 


power 
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An Expansion Method for Parabolic Partial Differential 
Equations 
J. W. Green? 


The aim of this paper is to adapt to certain parabolic partial differential equations an 


expansion method of solution developed by 8S 


to make possible a moderately compact 
the heat equation. Similarly, the 
general, but are assumed to be 
ean be reduced. The 
describes a definite 
bility of estimating the error. 


method of solution 


1. Statement of Problem 


In the region S, O<2<n, 0<t- , we seek a 


olution of the equation 
L(u) =uz,—Uy—gu f (1) 
ibject to the initial conditions 


u(r,0)=u(0,t)=u(r,t)=—0. (2) 
The functions f and g are continuous in S, together 
ith their first two partial derivatives. As a com- 
itibility condition, we require that f(0,0)=f(x,0)= 
is convenient to assume that g>0. This is not 
ill essential, since a substitution of the form 
e“r transforms (1) into 
ez—Vi— (9 +A)v mw. (1’) 
The new coefficient g+ can be made positive in 
s large a part of S as desired, and this is sufficient 
r our purposes. 
This method can be applied to more general equa- 
ons in which u, has a variable positive coefficient 
d the u, term is present. However, the equation 
ected is of sufficient generality to illustrate the 
iethod and results and yet not make the computa- 
fons unpleasantly complicated. In case g=constant, 
ie method reduces to the method of expanding in 
ize nfunctions. 


2. The Moment Method * 


We seek an approximation to u in the form **° 


i 
} 
| 
i 


> 


u(x,t) C,, +(t) sin ka, (3) 


ation of this paper was sponsored (in part) by the Office of Naval 
y of California and National Bureau of Standards 
lo, Un nuovo metodo per l’analisi esistenziale e quantitativa dei 
propagazioni. Annali della Scuola Normale Superiore di Pisa, [III], 
¥ (1047) 
ialerkin's method. See, I. S. Sokolnikoff, Mathematical theory of 
McGraw-Hill Book Co., Inc., New York, N. Y., 1446 
7 of the reference given in footnote 3, after the existence of the solu- 
proved, one can show that the system {sin &z} can be replaced 
table orthonormal systems. 
© 


Zand | means 


presentation, the 
most general to which the method is applicable, 
boundary values and initial conditions are 
in a canonical form to which others, 
not only 
procedure for approximating it 


Faedo * for hyperbolic equations. In order 
equations treated are not the 
but are the simplest nontrivial relatives of 
not the most 
if sufficiently smooth, 
existence of a solution, but 
remarks are made on the possi- 


shows the 
Some 


subject to 


S {L(u,) 


where 


4} sin jrdxr=0, 


Cc, (0) =0. 


n 
Equation (4) amounts to 


CO — SS Crabs (t) 


4 


—7?C,,;(t) a,(t)=0, 


where 


9 
= fg(x,J) sin jx sin krdz, 
T 


b, (0) 


9 


a,(t)=— ff (x,t) sin jrdz. 
T 


The functions C,.,(t) are thus determined by 
system of n simultaneous first-order linear differential 
equations. These equations are already in canonical 
form, that is, the derivatives are solved for in terms 
of the other quantities, and so they have a unique 
solution satisfying (5). Thus the functions u,(z,t) 
are determined in unique fashion 

We intend to show that the sequence u,(z,t) 
converges to a solution of (1). The method proceeds 
by showing that the wu, and certain of their deriva- 
tives are uniformly equicontinuous. Recall that the 
functions f of a family are uniformly equicontinu- 
ous in a region, provided that ¢« is given greater 
than 0, there is a 6>0 such that |f(p)—f(q)|<e 
for any f of the family and any two points p and q 
of the region whose distance does not exceed 4. 
Now, Ascoli’s theorem asserts that out of uni- 
formly bounded and equicontinuous sequence of 
functions a uniformly continuous subsequence can 
be selected. Therefore, by selecting subsequences 
of subsequences a sufficient number of times, we 
shall arrive at a subsequence of the u, which con- 
verges uniformly, together with sufficiently many 
of its derivatives to allow us to make the existence 
proof. This latter is effected by using a Green’s 
function and showing that the u,, for n>, come 





successively closer to satisfying an integral equation 
corresponding to (1), and thus that the limit function 
does satisfy such an equation. 

To prove equicontinuity of our sequences, we make 
use of a criterion of Tonelli’ as follows: A class of 
continuous functions p(z,f) in the region Sp:0 <2, 
0<t<7, is uniformly equicontinuous in this region 


if (a) for each ¥ and t, v(%.t) and v(z,t) are absolutely 


continuous functions of ¢ and z, respectively, (b) the | 


partial derivatives v, and v, are summable in Sr, 
(c) there exists a constant ¢ such that frdr<e and 
JS vidx<e for all members of the family. Further- 
more any accumulation function of a class satisfying 
(a), (b), and (c) also satisfies (a) and (b). 

Now the functions u, and those partial derivatives 
in which we are interested are easily seen to be 
continuous and to satisfy (a) and (b), and the 
burden of the proof of equicontinuity falls on obtain- 
ing inequalities of the form in (c). This is done by 
making use of devices similar to those frequently 
employed to show uniqueness, and to recall them, 
we shall first prove a uniqueness theorem. 


3. Uniqueness 


Suppose that two solutions of the system (1) (2) 
exist, all derivatives appearing in (1) being con- 
tinuous in S. Their difference v will satisfy (1), (2) 
with f replaced by zero. Thus L(v)=0 and so 


SL(v)vdr=S'{ va2—v,—gv}vdr=0. — (7) 
Now 


d 


Sv, edi S43 di v'da S v'dz, 


and 


S vervdi 


becomes 


S vdz, (7 
— J (v3 


If we put K(t) 


LK’ (t) -ge*)dzx <0, 


since g>0. Since by (2), K(0)=0, we see that 
K(t)<0. But K is by its nature nonnegative and 
K=0, which implies that v=0 and that the two 
solutions are identical. 

Theorem 0. There is at most function, con- 
tinuous together with all derivatives appearing in (1) 
which satisfies the system (1), (2). 


4. Certain Bounds 


one 


This section is devoted to obtaining the inequali- 


ties required to carry out the program outlined in 
section 3. 

4.1. Lemma. 
satisfying fer x>0 the 
where a>0. Then forz 


Let K(x) be a nonnegative function 
inequality K’(x)<ayK(z), 
0, yA(2) <yA(0)+(a@/2)z. 


? Tonelli, L'estremo assoluto degli integrali doppi, Annali della Scuola Normale 
di Pisa (1933). 





Suppose that for some positive 2%, , 
y K(0)+(a@/2)%. Then there 
VK (2,) =v K(0) + (a/2)x,, but y A(z) 
for 2,<2r<2. In particular, A(z) 


and so in that interval 


exists 2, su 
-yK(0) 
0 ing 


K'(r) 
= Sa 
\ K(x) 


Integrating from 2x, to 2, we have 


VK (2) — y K(x) < (a/2) (29 


or 


) 


VK (29) < yA (2,) + (@/2) (219 —2,) = VA(0) 


a contradiction. The lemma is proved. 
Corollary. Let K(x) be a nonnegative 
satisfying for r>0 the inequality 


K’ (x) < a(x) y K(2), 


where a (fr) is nonnegative and nondecreasing 
for x>0, 
VK(r) < VA(0)+ }2ra(z). 

Consider any %>0. For xr<2x, K’ (x) <a(a)\K 
By the preceeding lemma, , A(z) < y A(0)+ }als 
Setting z=z7, and noting that zx, is arbitrary, we ¢ 
the stated result. 

4.2. The bound® for fridr. Let each of 
equations (4) be multiplied by C,.,(é) and the res 
summed. There comes out the equation 


S {Uz 


U,—gu f udxr=0. 


Since 
S u,,udi S wdz and fuudr 


(10) may be written 


1 
a < utd 
2 dt J ’ 


S (e+ Sfudz< JS tud 


qu’ )dxr 


5 Swe 


Set /wi(z,t)\dr=K(xr). From (11) and Schwarm: 


inequality, 


LK’ (2) < S | fu dz <y S fdz-y S waz. 


maxy J f*(z,7r)dr= Fit). 


rst 
Then (12) gives A’ (t) <2F(t)yA(t). From the com 
lary of 4.1 we deduce that yA(t)<yA(0) +t/! 
Now K(0)= /w(z,0)dz=0, and so ¥ A(t) <t/(t 


* In section 4 we shall frequently omit the subscript n from t,. When 


u appears in section 4, it will mean u,. 
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Sw (2,t)dz<?F*(t), (14) 

s given by (13). 

e bound for fr2.dr. Let the equations (4) 
entiated with respect to t, multiplied by 
ind added. There results from this 
gu,—gu—f,}udr=0. (15) 


User — Ute 


. d . 
lr=— fu?,dz and fu,y,u,dt ba Sudz, 


r K(t)= fwadz, we get 


- I (u2, f,)u,dx (16) 


<S (\giul| +f, 


} qu? )dz S (guu | 


) iu, idx 


<V JS (geu Ff \)*dz-yK(t). 


max |g;(27,t) Gy (t 


rst 


’ 


max y S fi r.r)dz F(t). 
r<t 


Ie Vdr< S (Gi (t)w? ; 2G, t) uf; +f? dx 


G(t) fwdzrt+-2G, (Oy S wdzS fide+ Sfidz 


G? (t)O F(t) 4-24, (O)tF (t) F(t) + F(t) 
his last equation being a definition of //(t). 
‘rom (16) we have K’(t) <2H/(t)y A(t), and by 4.1, 


yA(t) <yA(0)+tH(t). (19) 


K(0) =f u? (2,0) = S(S5C,,.(0) sin kx)*dz 


U 
7 , _ 2ey re , 
= C.,.(0)) ~ (SS f(2,0) sin krdzy 


<S f%(x,0)dz 


y Bessel’s inequality. 
19) we conclude that ¥ A(t) < F(0)+tH(). 
Theorem 2. 
S 2, (z,t)dz < (F(0) +tH(t))?=A(t), 
is giwen by (13), 
H G3 (tHe F?(t) +26, ()tF (t) F(t) 4 


F?(t), (21) 


mi G, and F, are given by (17) and (18). 


FT? (t), | 


If in (4) one puts ¢ 


4.4. The hound of S @ de. 
written 


Equation 10 may be 


F —w—uu,—gu’ fu dz=0 


S wdz yy 


Sv JS wider: fujdr+ 


gu? +uu,+fusdr< f(\uu, 


ju da 
VS wdr-ffr*dz 


< VF? (t)A*(t) + VelF*(t) F(t) 
tF(t) (A(t) J F(t)). 


Theorem 3. 


Su (z,t)dz <tF() (A(t) + F(t) =Bb), (22) 
where A is give n by (20). 

4.5 The bound of SJ €@,.(z,t di We differentiate 
(4) twice with respect to ¢, multiply by C.,(t) and 
sum, to get 

£ Urrit 24; Jacl Fer Ue =O. (23) 


Une Gua 


As before, we write 


S User dle SJ ui,,dz, 


ar 
S Uh lz 5 it S u2, dz, 


S ui dz, 


| and (23) becomes, on setting A(t) 


K’(t) S 6 w,+gu}, }dz—Sf fie}u dz 


2 Uy rll T 


<J{ 29 uy Jit | 4 Su Ur \dx. 
It is clear that by the use of Schwarz’s inequality we 
can obtain an inequality of the form 
K’'(t) <Jid)vAK(o, (24) 

| where J(t) is an increasing and continuous function 
expressible in terms of G,, F, F\, and the similar 
functions G, and F, corresponding to gi fre. 

We estimate K(0) as follows. Differentiate (4) 
with respect to ¢ and put t=0, to obtain 


=S°{ Uzz,(2,0) 
= 


fF (x,0) }sin grdr. (25) 


S Uu(2,0) sin jrdzx gu,(x,0) 


Thus A(0) <F?(0), and from | 


for 7 Now w,,(z,0) is a linear combi- 
nation of sin jz, j=1, 2, ...n, and by (25) its 
first n Fourier coefficients are the same as those of 
Use:(2z,0) —gu,(z,0)—f,(z,0). Thus by Bessel’s in- 


1,2 n. 


equality, 


S 2, (2,0)dz < S { Use:(2,0) —gu,(2,0) —fi(z,0) }?dx. (26) 


0, there results 





— fu,(z,0) sin jrdzr= ff(x,0) sin jrdz, 
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Two integrations by parts gives 
— Sf tigr,(2,0)sin jz dr= ff,.(x,0)sin jx dz, 


and since u,;,(z,0) is a linear combination of sin jz, 
j=1, 2, . . ., mn, We see as above, 
Su, (z,0)dz< Sf f2,(2,0)dz. (27) 


By using Schwarz’s inequality, (27), and the bounds 
previously obtained, the right hand number of (25) 


can be dominated by a continuous function depend- | 


ing only on f and g and their derivatives. That is, 

K(0) is bounded uniformly in n, and from (24) it 
follows that the same is true of A(f). _ ' 

Theorem 4. There exists a continuous function C(t), 

de te rmined by fg and the ir de rivative 8. such that 

S08  (2dr<oC(t). (28) 

4.6. The bound for fr? dz 


From eq (15) we 
have on replacing /u,,u,dzr by —f 


S uidz, 


Lqutf,} ured. 


S uidz S User guy 
Because of the bounds previously obtained, it is 
clear that we have the theorem: 

Theorem 5. There exists a continuous function D(t), 
expressible in terms of f and q and their de rivatives 
such that 

S12 (2 ddze< Dit). (29) 

4.7. The bound for fw2_dr. If equations (4) be 

multiplied by 7? sin jr and summed, we have 


Sid2=S (upt+gu 


qu- f “dry S u2dz. 


-fhu,dr 


<vJ} ust 


Now from the previous bounds we can find a con- 
tinuous function /(t) which depends only on f,g and 
their derivatives, such that 
S {uet+gutf dr< E(t). (30) 
Thus we have the theorem: 
Theorem 6. There exists a continuous function E(t) 
independent of n such that 
Su. (ztdz< E(t). (31) 
4.8. An integral inequality. Consider the Green’s 
function G(z,£) relative to the equation y/’+dy=0 
y(0)=y(r). Specifically 


rir—eé 


G(z,é) 





and its uniformly convergent Fourier expsnsioy ; 


9 
aad Xs 

wiahe . 
a | 


I 


G(z,t) 2. sin jx sin jg 


For brevity, set 
u,t+gu+f=—*(u). 


then 
Urzr—O(u) = L(u) —f. 


Multiply (34) by G(z,£) and integrate with resp. 
to x to obtain 
S\ L(u) 


—ulét)— fb(u)G(z,b)dr I} G(x.8)d; 


Now 
S { L(u) 


F\G(x,t)dzr 
2 <= sin jt 
—=-S> —= S{L(u) 
he 2 ta 


f SIN jr 


2 = sin jé 
—=>> — L(u)—f } sin jrdz 
T n+l Jj 


because of (4). Now L(u)—f=u,,—®(u), and 


n 
Urs > C.(O/ sin jx; hence 
S Uz; sin jrdr=0 


oly oa 
n, giving 


for j 


Si L(u) 


. 2? = sin 7¢ 
4) G(z,8)dr >> se” SH (u) sin jriz 
: i 


F n+ 
By Schwarz’s inequality, 


S { L(u) 


f G(z,i)dz 


and 


>> (fu) sin jrdz) 


n+l 


< 2a (Su) sin jrdz)’< if Seu 


By (30), {#(u)dr< E(t), and so from (36) we ca 
conclude 


J c 2 l T= 
A - G(2,&\ < _* At 
I u) J Ae a dz =F \ 373 \ 9 I 


Referring back to (35) we have a theorem: 
Theore m 7. 
, QE (t 
—u,(&t)— fP(u,)G(z,8)dzr| < ; 
\ 3mn 





ee TSION » 
5. Existence Theorem 


now the discussion of equicontinuity and 
ion of Tonelli given at the end of section 2. 
of the inequalities contained in theorems I 
6, we see that the functions u,, u,,, and u,, 
formly equicontinuous. Thus by Ascoli’s 
we can pick out, by a suitable diagonal 
a subsequence {»,,} of {u,} such that 


Um (r,t) u(z,t), 


Umr(2,t)— v(z,t), 


Umi rt) w(2z,b), 

woach being uniform in each case. Because 
uniformity, one sees that v=0uw/dr and 
¥. The derivatives u,,, Uz; U, exist almost 

here and are summable. We consider 

S P(v_,) G(2,&) dz. 


v_,(&,t) - (38) 


In jrdz 


u»,, theorem (7) tells us that the expression 


Also (v,,)—>¢(u) 


und it nds to zero as Mm~o, 
so we have 


ormly, and G@ is bounded, 
tt)h— f(r, 


\G(x,é)dé 


-—ulé,t)— { (u)G(z,t)dx=0 


lhe integral §6(u)G(z,t)dz has a second deriva- 


ive with respect to z, namely, 


lence 


&(u) = —u,—gu—f. 
i. exists, is continuous, and satisfies 


—Uzz(é,t) + O(u) =—0, (39) 


hich is identical with (1). Since u obviously 
atisfies (2), it is a solution of the system (1), (2) 
mtheregion Sy. But 7'is arbitrary, and a solution u 
f (1), (2) exists in all S. 

We have shown that a subsequence of u, converges 
) the solution u, together with first order derivatives. 
tually the entire sequence u, so converges. If 
his were not the case, there would be a 
such that, for instance, lim Un 


n 


From the sequence u, 


subse- 
> u at 


juence U, 
| 


me point z,t. we ‘eat pick 


nother subsequence in the manner of this section, 
uch that this subsequence has for limit a solution v 
yf the problem. By the uniqueness theorem, u=p, 
nd so lim Un, (z,t)=u(z,t), a contradiction. 
e have 
Theorem 3. The 


sequence u,(z,t) has for uniform 


any Sy the unique function u(z,t) satisfying | 


1 (2). Furthermore, u,,-u, and Uyz;>u,; 


ly in Sr. 


6. Convergence in Mean of u,,, 


The solution u(z,t) of the system (1), has a 


Fourier series 


u(x,t) 


>5A,(t) sin jz, 
l 





the | 


Thus | 


which may be differentiated twice with respect to z, 
and the resulting series converges in the mean to u,,. 
Set 


n 
R ' 
25 i(t) sin jr. 


(40) 


If (4) are multiplied by and 


summed, we get 


*(C'n s(t) —Aj(0) 


. 0° 
S { L(u,)—f} < (un 


or , dx =) 


(41) 


Since L(u)=f, we may write (41) as 


S {L(u,—v,)—L(u 


Recalling that 
— u)- f, 


L(u)=u,z2—U, gu—f-+ Bos 


we write (42) in the form 


Shea 
e -P(U,— UV») 


and Schwarz’s inequality gives 


+ f—L(u 5(U,—v,) dz, 
J 


es 
V Fira ee Un)¢ da 


<vJS{—P(u, Liu—v,)}*dz. (43) 
Now, —®(u,—?,) +f (Un —Vn) -—G(Uy,— 0») 0 uni- 
formly, and L(v—v,) tends to zero in the mean; thus 


the right hand me aie of (43) tends to zero and so 


S ses _ ra) } da >), 


From this and an appropriate triangular inequality 
it follows that 


Theorem 9. The second derivatives Unzz of the 
approximating | functions tend in the mean to Uz. 

It is too much to expect that u,,,—u,, uniformly, 
or even pointwise. For example, u,,,(0,t) =0, w hile 
u,,(0,t)=f(0,t), and this latter does not need to 
vanish identically. 


7. Estimating the Error 


From theorem 9 it is easy to prove that L(u,) >L(u) 
in the mean. Set L(u,)=f,. Then 


im. 


L(u,—u) = L(u,) — L(u) =f, —f 0. 








Using the method of 4.2, we prove that 


S (u,—u)*dz <P F(t), 
where 
(44) 


F(t) =maxy J { f(z,r) —f, (2,7) }*dz 
rst 


Theorem 10. The mean square error satisfies the 
inequality 

S (u,—u)*dze <0F? (45) 
where F%(t) is given by (44), and Frit 

Theorem 10 unfortunately gives information only 
about the mean square of the error instead of about 
the error at a point, which would be desirable. This 
latter will be mentioned below. Apart from this, 
however, theorem 10 gives a rather curious sort of 
information about the m.s.e. It does not give an 
a priori estimate of this error; that is, we cannot say 
that after n steps the error will have a certain bound, 
because we do not know how rapidly /, tends to 
zero. On the other hand, we do seem to know more 
than the mere fact that the m.s.e. tends to zero. 
We know that if we take n sufficiently large, the 
m.s.e. will be small and we shall know it to be small. 
What we do not know is when we shall know it! 

In section 5 it was mentioned that u,, is a sum- 
mable function. Thus from L(u)=f, it follows that 
User 18 Summable. From the equation L(u,—u) 

f,—f, we obtain, using the methods of sections 4.3 


and 4.4, the inequality analogous to (22), 


S { (us 


where FY and A® are given by (13) and (20), except 
that f,—/f is used insted of f in their definitions. 
Now we are unable to say much about the factor A* 
since it contains terms of the form /{f,,—f,}*dz; 
that is, depends ON Use, ANA Unsere. However, it 
appears quite likely in many cases that A‘ will 
remain bounded asn—. In any case the inequality 


ds no, 


u), }*dz <tF* (t)(AS(t)+F*(t Ge®, 


S { (u,—u), }°dz < G* (t) 


is a valid one, and in a given computation G* 
actually be computed. Setting A,(z,t)=u, 
have from Schwarz’s inequality 


be 3 ~ 
A, (2,t) | Anz(s,t) ds Sa | A,2(s,t) }*ds (4 
* ) . 


<V tJ { Ansls,t) Pds<V 2G, (0. 


from 7 instead of from zero, 


A, (2,t) <V(4r—2) GF (t). 


From (46) and (47) it follows that 


9 


A, (z,t) Sy; a(x—x)G* (0). 
T 
Theore m il. 


u,(2r,t)—u(z,t) 


whe re 


Gi()=tFS(){ AXO+ F320 


F*(t)=max y J { f.—f }*dz, 


r<t 





A*()—F3(0)—tH* (0, 


H(t) G(OlC { FF (0) }?+- 2G, (OFS (O+ FE tt 


In 


Fs. 


(O)=max y J{ far(z,7)—f(z,7) }2dz. 


rat 
For any stage n of the solution, the function @*! 
| can be explicitly calculated. If the A} do remain 
bounded, then G30, and we will get a_ useful 
estimate of the error. 


Los Ancretes, April 24, 1952. 













In an analogous fashion we obtain by integrating 
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Characteristics of Internal Solitary Waves’ 


Garbis H. Keulegan 


This is an application of the method of approximations initiated by Boussinesq to the 
disturbances of the interface points for waves of po form and the internal solitary 


wave 


wave is established 
and on laver thickness is determined. 


1. Introduction 


nvestigations on the model laws of density 
is, a project now being carried out at the 
ial Hydraulics Laboratory, need arose for the 
ration of the genesis and the damping of 
il waves at the interface between layers of 
ind saline waters. Among the many possible 
of such disturbances, one may consider for 
ake of simplicity the behavior of a single 
escence and of nearly sinusoidal progressive 


has been our purpose to deal with these two 

of internal wave motions experimentally. 

experimental studies of solitary waves have 

completed. As a natural guide in the study 

data, resort has been made to a theoretical 

The present paper gives the basis and 

result of this analysis. Relations are here 

obtained giving the dependence of the velocity of 

wave propagation on wave heights, the form of the 

olitary wave, and the expressions for the velocity 

vector on the upper and the lower layers. Con- 

sideration of the experimental data, however, is 
reserved for a future occasion. 

After putting forth the basic conditions for the 
analysis, the question of internal waves of infinitesi- 
mal height and negligible interfacial surface curva- 
ture is taken up in the first approximate solution. 
In the second approximate solution the character- 
sties of solitary waves are revealed. In general, 
there are seen to be some similarities between 
ordinary surface solitary waves and internal solitary 
waves. No attempt is made to extend the approxi- 
mations to a third-order analysis as no special de- 
mand is made by experimental evidence for this. 


2. Mathematical Formulation 


A layer of lighter liquid of thickness Hl’ and of 
y p’ rests on a layer of a denser liquid of thick- 
ness // and of density p. The upper liquid at its 
iree surface is exposed to air and the lower liquid 
rests on a rigid horizontal bed. Both liquids are 
initially at rest. The displacement of the interface 
with respect to its initial undisturbed position is 
denoted by h; the displacement of the free surface 
with respect to the level of the undisturbed free 
surface is denoted by A’ (see fig. 1). Taking the 


—_ 
t 


densit 


aration of this paper was sponsored in part by the Office of Naval 


Sear 


The system considered is a layer of liquit 
liquids of the layers being initially at rest and of constant total depth. 
rhe dependence of wave velocity on wave height, on density differences 





on another layer of greater density, the 
The form of the 


origin, O, of rectangular axes at the undisturbed 
interface, the axis of z is drawn horizontally. The 
axis of 2 is drawn vertically and the positive branch 
points upwards. The velocity components along 
these axes are denoted by u, w, and wu’, w’, the 
primed symbols referring to particle velocities in the 
upper liquid. 

It will be assumed that the disturbanees are pro- 
duced in liquids initially at rest, so that the conse- 
quent flow is irrotational and admits the velocity 
potentials @ and ¢’. The vorticity that is naturally 
present at the interface will be ignored. Since the 
type of disturbance visualized is translational, that 
is, the particle velocities in vertical planes nermal 
to the direction of wave motion are nearly constant 
for each layer, it is appropriate to introduce the 
expressions of the velocity potentials 


: oO 
v7) (cos z 5 yt ‘s. ) (1) 


for the lower layer, and 


+) ¢ ( sin 2 


for the upper layer. Here, o, 4%, 
functions of x and ¢ alone. 

Altogether there are six unknowns, h, h’, $0, 9, 
4, and 6, which are determined on the basis of the 
kinematic and dynamic conditions at the various 
boundaries. The pressure is atmospheric at the 
free surface. At the interface the pressure is con- 
tinuous. A particle in the free surface remains in 
this surface. The particles in the upper surface of 
the lower liquid at the interface remain in this 
surface. Finally, the normal particle velocities 
vanish at the horizontal bottom. 


< 


H 


¢’ ( COS 2 


do, and 6, 

















° 
H P= P'+AP,$, 
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Ficure 1. Notation diagram. 
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The mathematical formulation of the dynamic 
boundary conditions requires that expressions be 
given for the pressures at all points in the two 
layers 
reduced to Zero, 


P _. (u?-+ 1° oH’ 
p ol ; 2 i. 


for the lower liquid, and 
(4) 


Translating now 


for the upper liquid. the state- 
conditions 


ments of the boundary into 
mathematical relations, neglecting w* with respect 
to u*?, and w” with respect to vu”, they are 


above 


09". bans u” 
an 3 


oh 


ol 


and 
(10) 


| Hence, 


Assuming that the atmoshperic pressure 1s | 


¢ oo—( Hz + ~) oe, 


and therefore 


Op _(H+2) Od 


Oz Or 
Similarly, from eq 2, 


* 0’, 


dot? -3 or’ 


or 
hence 
0g’ 086, 


_ , bo 
Oz Or 


Or" 


Taking the kinematic boundary condition of 
| upper surface, eq 6, and neglecting the produ 


3. Internal Waves of Infinitesimal Height | 


The functions @ and @’ may be expressed in 
terms of @ and @»’, respectively, after noting the 
kinematic boundary conditions at the bottom and 
at the free surface. From eq 1, keeping only the 
first three terms, 


_ 06) 2° O'bo 


Oe 3 Oz" (11) 


P= P07 


and therefore 


Oo 08, _ Odo 
Oz Or ' OF 


The latter in view of eq 10 requires that 


OO) Od» 
Or ~H or? 
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| 
| 
| 
| 
| 


| 


term which is a small quantity, 


oh’ 


06, 0d 
= , ] hh’ 
Or (H sh Or ol 
And as h’ 


may be neglected in comparison with H 


06> = oh 4 H’ O"do 
Or ol Or’ 
hence 


and therefore 


0d’ 
or 


) O° go, 


- ae 
Se FH —2) S33 


For the points at the interface, that is z=A, o1 


may write from eq 14 and 15, neglecting the terns 


involving z and its square, 


Og Ody 
ot of 
and 
Og ae O° >» 
oz H oz 


From eq 20 and 21 
Oo’ _ Ogp 
ol ol 
and 


og’ _ ° h’ , O’do 
= | +H oz? 


| These imply that wave height A is small in compari 


with H or H’ and also that the curvature of the wa’ 
surface is small. 
Taking next the boundary conditions, eq 5, 7,° 


and 9, neglecting in them terms consisting of square 





rhis 


al 


and 


parisol 
e wal 


square 


oducts, and using the relations from eq 22 to 


usive, 


Od; 
OL 


gh’, 


p’ ( Odo of ) p ( OFe_ op ), 


ol 


H 


O* do 


—a Oz 


ot 


Od» Oh 


or? Ot 


_ “4H? 29%, 


(29 
a Oz* 


inging the terms in eq 27 and then differ- 
ng twice with respect to t, we have 


0°¢, p O° 


oo p—p’ , oh 


oe! p’ de p »* of 


ntiating eq 27 twice with respect to gz, 


0'¢, p Odo p c. = he 


— (31) 


1022p’ O02? p da 


tiplying eq 29 by g, introducing the value of h’ 
eq 26, and differentiating ‘the result with 


respect tot 


—gH 


ozr°0t 


» Go ___—O*ho 
gH’ + 291 or 


[his reduces, after eliminating ¢, by means of eq 


and 31, to 


O°, 1 P °° 
- xo gH ‘Ox 


do ,(p—p’) oh 
+0°H ; 
“Ol we p Ox 


p Od p—p’ Oh 
p’ ot° + p’ g ot* 


Differentiating the latter with respect to ¢, and then 


eliminating A by means 


0'¢ = do 


g(H+H’) Or 


Putting p—p’=Ap, 


of eq 28, there results, finally, 


) . F o* ) 
+9°H’H 2&—? - 0, (32) 
p 


which is the differential equation of wave motion in 
a liquid system consisting of two strata of different 
densities for small disturbances. 


wig H'H 


w= 
i 


@) Tw 


gH+H’), 


ve equation, eq 32, may be rewritten 


°* o? o’ . & 
(sa-“t on :) (Sa-“ 3 5z3) H=0. 


Two types of waves are possible. In one the waves 
have the velocity of propagation +); in the other, 
the velocity +... The plus sign refers to waves 
moving in the direction of z negative and the minus 


sign, in the opposite direction. Furthermore, these 


waves, of infinitesimal wave height and of negligible 


| surface curvature, travel without distortion of form. 


The discrete values of the velocities may be shown 


| from eq 34, neglecting the higher powers of the 


density differences, to be 


H’'H be) 


Wt rite +-H] (1 HH? 


and 


gH’ 7 pe a H Ap 37) 
H- +H’)? p ae 
Of these, the second refers to disturbances of the 
interface, that is, to internal waves. The first refers 
to ordinary waves; in this expression one notices the 


| effect of the nonhomogeneity of fluid on the velocity 


of propagation of the ordinary waves. Nonhomo- 


| geneity reduces the value of the velocity of propaga- 


tion. For a given relative density difference the 
reduction is greatest when the layers are of equal 
depth. For the purposes of the present work, the 


secondary effects of the relative density differences 


will be ignored. Accordingly the velocity of propa- 


| gation of internal waves of infinitesimal height, 


replacing subscript 2 by 0, will be given as 


, gH’'H Ap 32 
o~H+H’ p- —— 

The next question to be considered is in regard to 
the particle velocities in the two layers. One com- 
mences the analysis with the lower layer. As was 
noted above, the internal waves of infinitesimal 
height progress without change of form. This fact is 
equivalent to the equality of the operators 


re) 
—=Fo- (39) 
ot , 

The negative sign is chosen for waves moving in the 
direction of z positive and the positive sign for waves 
moving in the direction of z negative. For the sub- 


| sequent analysis it will be assumed that the internal 


waves move in the direction of x positive, and thus 


d d 
ue Or (40) 


One now naturally selects that kinematic condition 
of the interface stating that a particle of the lower 
liquid once on the interface remains on the interface. 
This condition is given by eq 28, which now may be 
written as 


oh 


HS Ga-e 





and, as in terms of the potential, the particle velocity 
Lo is 
Od 
mo = J (41) 
Or 


we have 
OUy Oh 
TH] Ww . 
Or Or 
Integrating the latter equation with respect to z and 
observing that vw vanishes when A vanishes, then 


on (42) 
H 

Accordingly, particle velocities in the lower layer are 
proportional to the velocity of wave propagation and 
vary with the wave-element height of the internal 
wave. When these heights are positive, that is when 
the wave elements are elevated, particle motion in the 
lower layer is in the direction of wave motion. 

In terms of the velocity potential the particle 
velocity in the upper layer is 


Abo. 
Or 


To obtain the values of u, one now considers the 
other kinematic boundary condition of the interface. 
The reference is to eq 29, which when treated in an 


analogous manner yields 
H’u,+Hig=oh' (44) 


Now, when in eq 26 the partial differentiation with 
respect to ¢ is replaced by the partial differentiation 
with respect to z, the result is 
wy gh’. (45) 
Eliminating A between.eq 44 and 45, and introduc- 
ing the value of w from eq 38, one obtains the 
relation 
Ap H 


> H+’ (46) 


H'u, (1— )\+Huo=0, 


and this connects the particle velocities in the two 

layers. As the term containing the relative density 

difference is small, it will be neglected, and therefore 
u, H 


a (47) 


Thus the motion of the particles in the upper layer | 


is oppositely directed to the motion of the particles 
in the lower layer. 

It will be instructive to show to what extent the 
internal waves do affect the free surface. This will 
be better understood if the ratio h’/A is investigated. 
The combinations of the expressions in eq 42 and 
45 gives 


h’ wow uu, 


h gH uy 


> 


which, in view of eq 38 and 39, reduces to 


h’ _ Ap H 
h »p H+H’ 


The interpretation is that the displacement of ¢/ 
surface is directed oppositely to that of the interfa 
Furthermore, the disturbances of the free surfg 
are very much reduced in comparison wit! 
displacements of the internal waves. 

Finally, the question of the energy of in erhal 
waves of infinitesimal heights may be considere 
The energy of waves is in one part potential and 
another part kinetic. Evaluating the potent 
energy with respect to the undisturbed configuratio, 
of the two layers, the appropriate expression is 


— { ~_ U 


, he [*.. »’ 
E p=9 a h?dx- qd f 
v7 } 4 
where \ is the effective wave length, that is ¢| 
length along which A? is measurable. Expressing 
h’ in terms of A through eq 48, 


» Ap Ap H 9 weer 
Ep I> [1+ Gaze) If, hedz. 45) 


The kinetic part of the energy in terms of t 
particle velocities in the two layers is 


a.) 


E, a | uy H'de +8 | wH dx. 


Expressing uv in terms of uv with the use of eq 4 
expressing U» in terms of A and w with the use of 
42, and introducing the more exact values of w fro 
eq 37, the final value of kinetic energy is 


- Ap Ap H wit 
E, g 5 [ i+ > (7 i w) \J, h?de. 


Comparison of the expressions for the two energie 
shows that, in internal waves of infinitesimal wa 
height when the upper surface is free, the potentia 
and kinetic energies are of like value. When t! 
relative density-difference term is ignored the energ 
of internal waves is simply 


*» 
E gdp | h?dz. 
0 


The form of the internal waves together with 
difference of the density of the two layers is sufficien! 
for the evaluation of the energy. 

Some of these properties of internal waves of larg 
wave length and of infinitesimal wave height a 
well known. A short derivation of these results, {« 
example, by the method of Lord Rayleigh is give 


by Thorade [1]. 
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= osity. 


4. The Internal Solitary Wave 


neral when the wave-element heights of the 


1] waves are finite, a deformation of the wave | 


occurs during the travel of the waves even 
he two liquids are assumed to be perfect and 
The extent of the deformation depends on 
eht A and on the surface curvature 0°h/dz°. 


responding mode of the deformation in the | 


ry translation waves, as first discussed by 
nesq, is well known [2, 3]. A similar theory 
worked out also in the cases of internal waves, 
iis theory could form a basis for the study of 
al waves which travel without deformation. <A 
| procedure is to assume the existence of internal 
vy waves, positive or negative, which travel 
listances without deformation in the absence 
These intumescences will be referred 
nternal solitary waves, and the analysis to be 
ed below will determine the form of the waves. 
inarily in establishing the character of the 
ssive waves of the permanent type resort is 
to an artifice [4]. By superposing on the flow 


ent of the magnitude equaling the velocity of | 


gation of the wave and moving in the direction 


site to that of the wave, the system is reduced 


steady state. 


Another method, and this will be 
here, is that the time differentiations will be 
into space differentiations for the wave 


eling in the direction of z positive using 


fy 
ol pro 


The 


(40a) 


ere w is the velocity of progression of the wave of | 
rmanent type. 


In the present case w is the velocity 
eression of the solitary wave. 

expression for the vector potential in the 
layer to the second approximation, is, from 


0% 2° Od 2° 0°0 na 
12 — -- —-- . (52) 
= Po os 20s" 6 OF 


ating the last term on the right-hand side, 
¢ 


t is the term of smallest value, from the first 


oxiumative value, 


06 
Or 


O° do 
—H or ) 


tor potential now is 
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In view of the bottom kinematic boundary condition, 
eq 10, 
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which than is the second approximative value of 
06,/Oor. Thus, the resulting form of the vector 
potential is 
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This is the expression for the velocity potential 
in the lower layer to the second approximation 
involving the single unknown function ¢@. When 
the height of the internal solitary wave is smaller 
than the depth /7 of the lower layer, the expression 
for the vertical component of the particle velocities 
at the interface may be simplified by neglecting 
2 in eq 55. Mathematically the basis of the ap- 
proximation ts 
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Now, the last two terms on the right-hand side being 
the smallest in value, substitutions may be made 
for them from the first approximative values, 
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Introducing these in eq 55 and putting h, the 
vertical component of the particle velocity at the 
| interface, sign reversed, is 
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Similarly the expression for the vector potential 

in the upper liquid, from eq 1, is 
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The last term on the right-hand side being the'smallest { and 

term, the value of 6 may be taken from the first 

approximation, | a O60, H” o'¢, 2°H’ dt, 

; ; s = °d 2 Oz 2 d2* 

0% _oh 4H’ O"do. (19a) 

Or ot Oz | When the latter expression is applied to the reg 

a — rec i i close to the interface, A? being smaller than H’ 

It is permissible at this time to affect a small modi- | (oi involving 2? may be neglected. Thus, 

fication in the analysis, and this consists of neglecting . 

the term dh’/dt. The assumption implies that | ae! a ae a, , H” ate) 
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Introducing this is eq 55, Consider next the two kinematic boundary con- 
er — P | ditions for the interface given by eq 8 and 9. Thes 
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which is the second approximative value of 06//dz. 
Substituting in eq 61, 
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eq 73 the value of w’ from eq 69, the 
boundary conditions reduce to 
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The dynamic condition for the interface given 


eq 7 may be written in the form 
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The last two terms on the right-hand side are 
mall quantities and these may be evaluated by the 
first approximative values of the particle veloci- 
ties u’ = —woh/H’ and u=wA/H. 
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Expressing the value of the density difference 
Ap in terms of #3, using eq 38, and dividing 
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Multiplying by 0//dz, integrating once with respect 
to’x, and putting the constant of integration equal 
tofzero since all disturbances vanish at infinity, we 
obtain 


(S-1)i 


Let the crest height of the wave be fy. 
maximum, O0//Ox vanishes at this point. 


w* , (H'-H 
(=~! ) hi—( ) hi 


WH 
H’—-H , 
HH “ 


H’—H 


Oh ) 
HH 


Or 


0. (S86) 


)n— + WH ( 


As this is a 
Hence 


Ww Wo T . i 

\ H'H 
This is the law of velocity of propagation of an inter- 
nal solitary wave. The magnitude of wo, the velocity 
off{propagation when the wave is exceedingly small, 
Substituting in eq 86 from eq 88, 
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is given in eq 38. 
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and the solution of the equation is 
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This gives the form of internal solitary waves. 

Examination of eq 86 reveals that the relation o 
the thickness of the two layers has an importan 
bearing on the formation of internal solitary wayes 
When the depth of the upper layer is greater tha 
the depth of the lower layer, H’>H, the intern; 
solitary wave is of positive type, that is, A is positiy 
When the depth of the lower layer is greater {hg 
the depth of the upper layer, the internal solitg; 
wave is of negative height, that is, A is negatiy 
everywhere. When the two layers are of the sam, 
depth and the difference of the densities is very sma] 
the formation of a solitary wave is excluded by this 
analysis. 

The consideration of a fairly large number , 
experiments tending to verify the above theoretic, 
results, first in regard to the dependence of velocit 
of propagation on wave height, and second in regard 
to the form of internal solitary waves, will be resery: 
for another occasion. 
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Becker Value of Manila Rope by Photoelectric 
Reflectometry 


Sanford B. Newman, Harry K. Hammond, III, and Helen F. Riddell 


Becker value, or reflectance under specified conditions, has been used for more tl an 20 


vears as an index of fiber quality in the purchase of manila cordage, 
a visual reflectometer. 


for determining Becker value requires 


The accepted method 
Previous attempts to usé 


photoelectric reflectometers were unsuccessful because the instruments lacked sensitivity, 


attributable in part to the low transmittance of the chosen filter combination 


Properly 


selected filters used with a sensitive reflectometer permit photoelectric determinations to 


be made with greater speed and precision than visual determinations 


The results obtained 


photoelectrically agree well with the averages of visual determinations 


1. Introduction 


nila fiber (abacdé) is one of the important raw 
als or cordage. Rope made from this fiber 
msumed at an annual rate of over 50,000 tons 
9' {1)2 The Federal specification for manila 
2] contains requirements for strength, weight, 
‘ther physical characteristics. Opimion in the 
try, however, holds that serviceability is related 
or, the dark grades of fiber being less serviceable 
than the lighter [3]. Becker and Appel, at the request 
f the Cordage Institute and the Federal Specifica- 
soard, investigated the color characteristics of 

la rope fiber and in 1933 reported a quantitative 

al method of evaluation [3, 4]. One result of 
investigation was the establishment of the 

cer value” of manila rope, defined as 100 times 
ratio of the reflectance of a suitably prepared 

le of fiber to that of a standard magnesium oxide 
for monochromatic light of approximately 500 
mu. Fiber and standard were equally illuminated at 
with CLE (Commission Internationale de |’ Eclair- 

ge) standard source A and viewed perpendicularly. 
\feasurements were made with a Martens (visual) 
photometer with a Wratten filter No. 75 over the 

eDLEe¢ eC. 

Despite the lack of data on correlation of color 
vith serviceability, measurements of Becker value 
continue to be made and appear to have facilitated 
the purchase of satisfactory manila cordage. With 
the passage of time, considerable confidence has been 
placed in the efficacy of Becker value measurements, 
md they will probably continue to be used in 
cordage specifications for some time. 

Since the method was first described, Becker value 
has been determined by means of a reflectometer 
itilizing a Martens (visual) photometer. The pro- 
edure requires four pairs of observations, or settings, 
with the visual photometer, averaged to obtain one 

ination. As with all visual photometers, the 
‘e dependent on the skill and experience of 


Oo ct 


samp 


vel 


ent data have not been released : 
brackets indicate the literature references at the end of this paper- 


the observer. Furthermore, the Martens photom- 
eter scale is graduated in degrees so that Becker 
values cannot be obtained directly from readings of 
the photometer. A precise direct-reading photom- 
eter would be more desirable. Thus there arose 
the desire to replace the visual apparatus for deter- 
mining Becker value with photoelectric equipment. 
Therefore, a photoelectric reflectometer * was adapted 
to meet the requirements for Becker value measure- 
ment. 


2. Instrumentation 


The optical arrangement of the photoelectric 
reflectometer chosen for use in determining Becker 
value is shown schematically in figure 1. The filter 
system first tried consisted of (1) a Barnes luminosity 
filter [5] to correct the response of the photocell to 
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Figure 1. Schematic diagram of photoelectric reflectometer 


8, Cell containing cut fiber 
lenses; L, light source 


F, filters; P, photocell; M, mirrors; C, condenser 


! This is the photometric unit with 45°0° exposure head designed by 


Hunter 
(Gardner Laboratory, Bethesda, Md.] 
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FIGure : 


ectometer 


in the photoe lectric re 


that of the CLE standard observer, and (2) a Wratten 
filter No. 75 to obtain essentially monochromatic 


light at 500 my, as required for the measurement of 


Becker value. Although nominally duplicating the 
spectral response of the visual method, this combina- 
tion produced reflectometer readings considerably 
higher than those obtained visually. Substitution 
of a Hunter (bluer) luminosity filter [6] for the 
Barnes filter produced somewhat better agreement 
between the visual and photoelectric values. Inves- 
tigation of the spectral characteristics of the source- 
filter-photocell combination revealed that infrared 
radiation was responsible for the high reflectance 
values obtained with the photoelectric reflectometer. 
Figure 2 shows the spectral response of a barrier- 


layer type photocell and the spectral transmittance | 
The Barnes and Hunter | 


of the Wratten filter 75. 
luminosity filters also transmitted some infrared 
radiation. To eliminate the infrared component 
and to obtain maximum sensitivity at 500 my, the 
luminosity filter was removed and a Corning glass 
4303 was used to absorb the red and infrared radia- 
tion. This glass in combination with the Wratten 
filter 75 gave quite satisfactory results, as will be 
shown later. 

The specimen holder ordinarily used with the 
Martens photometer provides too small an area of 
sample to permit the use of the entire illuminating 
beam of the photoelectric reflectometer. The speci- 
men cell shown in figure 3 was constructed to permit 
large area illumination of the sample in the photo- 
electric reflectometer and still permit measurement 
in the visual instrument. The cell is constructed of 
methyl methacrylate. This material is easily worked 
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Plastic specimen cell used in photoelectric dete 
nations. Ol 


Ficure 3 


A, Body of cell; B, back of cell. The back is held in place by the t 


rubber bands which are looped over the protruding studs in A and contact th 
back on the rectangular raised portion, 

yl 
with ordinary laboratory tools, and the plastic pro- par’ 
vides a good, durable substitute for the “‘water white not 
cover glass, which often is not readily available. Tly Fro 
plastic cover is softer than glass and will abrade wit! tu 
use, but replacement is easily accomplished. A, perl 


electrostatic charge, induced by cleaning or filling th: rope 
cell, can be removed with a gas flame or high relativ: 
humidity, as is done with conventional cells. Dif. 
ferential orientation of the fiber caused by a charge grea 
cell is readily detected photoelectrically by observing 
the instrument galvanometer while rotating the ce! A 
in its own plane. Specimen heating by the ligh 
source is much reduced with the photoelectric instrv- 
ment as the light source draws about 20 watts as 
compared to 1,000 watts for the source in the visual 
instrument. 

As with nearly all reflectometry, the primar 
reflectance standard for Becker value measurements 
is a freshly prepared layer of magnesium oxid 
Working standards of porcelain enamel, however 
are much easier to use. For greatest accuracy wit! 
visual or photoelectric reflectometers, a working 
standard should have spectral characteristics sim- 


cius 
Bech 
Irom 
that 
sary 
ene 
the < 
show 
same 


lar to those of the specimens being measured. Becker Ph 
reported the spectral reflectance of fiber from 25 sam- are | 
ples of 1-in. manila rope from 13 manufacturers that 
and the spectral reflectance of the Philippine Island lew t 
Government standards [4]. The reflectance at 500m, In 
ranged from 20 to 60 percent. The spectral reflee- elect 

Visua 


tances increased nearly linearly from 450 to 550mz « 
the rate of approximately 0.15 percent per millimi- 
cron. A pale-yellow enamel working standard wit! 
a factor of 55.0, as related to 100 for MgO, was use 
in most of this work. 
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3. Observations and Results 


es of rope that had been submitted for test for 
mance with the Federal specification [2]. The 
vas prepared for test as described in the stand- 
ocedure, except that a rotary knife was used 
itting the fiber into the short lengths required. 
cnife was driven slowly, and no detectable 
ing of the fiber occurred. Reflectance meas- 
ents were made visually and photoelectrically 
it disturbing the contents of the cells. In fig- 
the values obtained visually are plotted against 
obtained photoelectrically. The relationship 


ker values were determined on several dozen | 


en the two sets of values is practically linear. 
quently, a straight line was fitted (broken line 

1), using the method of least squares, and 
dering the precision of the photoelectric values 
preciably better than that of the Martens values. 
quation for this line is 


+0.972, 


y=1.21 
r is the photoelectric value and y the Martens 


tatistical analysis of the data shows that the de- 
ire from a zero-intercept unit-slope relation is 
attributable entirely to random fluctuation 
a practical viewpoint, however, these depar- 
tures are unimportant. The minimum Becker values 
permitted by the Federal specification [2] are 40 for 
rope over 2 in. in circumference, and 43 for rope 2 in. 
reumference and under. In the Becker value 
range 35 to 55, the range of practical importance, the 
greatest discrepancy between the two methods is 
only 0.3 

\ statistical analysis has been made of 25 determi- 
nations of Becker value by one operator using one 
Martens photometer. From this analysis the con- 
clusion was drawn that a single determination of 
Becker value by the visual method is unlikely to be 
error by more than one unit. Comparison of data 
from different laboratories [7], however, indicates 
that as many as four determinations may be neces- 
sary for this precision. Observers of unequal experi- 
ence show considerable variation of Becker value on 
the same specimen, and even experienced observers 
show a range of about 1.5 for 25 observations on the 
same sample. 

Photoelectric determinations, on the other hand, 
are not influenced by the operator. It is estimated 
that instrument variations amount to as much as a 
few tenths of unit of Becker value. 

In conclusion, it may be stated that the photo- 
electric method is considerably more precise than the 

al method and that the two methods show a | 





€ 


Becker values 


values 


Ficure 4 Photoelectris plotted against vis sal 


Broken line represents the least-square line, 


systematic variation of not more than 0.3 unit over 
the extreme range of interest. 


The authors acknowledge the assistance of John 
Mandel in analyzing the data. 
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Dielectric Relaxation in a Styrene-Acrylonitrile Co- 
polymer During and After its Polymerization” 


Paul Ehrlich* and Nicholas J. De Lollis 


Dielectric relaxation in plasticized copolymers of styrene and acrylonitrile is studied 
during their formation as well as in the fully polymerized copolymers by measurements of 


dielectric constant and loss factor at several frequencies from 100 c/s to 100 ke/s 


In the 


polymerizing mixture, there occur, at an early stage of the reaction, sigmoid decreases in the 
dielectric constant at each frequency, accompanied by maxima in the loss factor, when these 


variables are plotted as a function of reaction time 
resulting from the relaxation of the nitrile groups. 


These changes are interpreted as 
The electrical properties of the fully 


polymerized copolymer, as it goes through its glass transition, are in semiquantitative agree- 


ment with those of the polymerizing mixture at the stage of the reaction referred to, 


dem- 


onstrating the occurrence of similar phenomena in each case. 


1. Introduction 


0 
roperties throughout a polymerization reaction have 
en reported. Lazarev and Raschketaev |1]‘* have 
oted a dependence of d-c conductivity, dielectric 
onstant, and dielectric loss on the degree of polymer- 
ation of the product resulting from the heat treat- 
ent of drying oils. Majury and Melville [2] have 
drawn certain conclusions about the nonstationary 
in vinyl-type polymerizations from the observa- 
tion of very small and very rapidly occurring changes 
n dielectric constant observed in the early stages of 
the reaction. This paper reports on a direct study 
of changes in dielectric constant and dielectric loss 
occurring throughout the polymerization of a polar 
casting resin. The measurements were carried out 
while the reaction mixture was polymerizing and 
ross-linking from a viscous solution containing 
mostly monomer to the fully cured polymer, which 
was a hard solid. 
According to the Debye concept of dipolar relaxa- 
}} the contribution of dipole relaxation to the 
, Is given by 


; 


tion 

lielectric constant, ¢’(w)—e, 
; 4 > , 
€ (w)—e, =(€,—e,,.)/1+a*r’, 

and the loss factor, e’’, resulting from this process is 

riven DV 


s7 , > 
€ '(w)=(e€,—€,)wr/1-+a*r*. 


Che loss factor thus goes through a maximum at a 
irequency fmax, Such that 


r=1/22f max: (1) 


Furthermore, the Debye theory predicts for a dilute 

suspension of dipolar spheres in a nonpolar solvent 
. 1 by the Department of the Army, Office of the Chief of Ordnance 
et lin part before the American Physical Society, Durham, N. C., 
i 8, 1953 

., Present address, Department of Chemistry, Harvard University, Cambridge, 


"f brackets indicate the literature references at the end of this paper. 
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that 
4xrna’ kT, 


where « and e¢, are the static and infinite fre- 
quency values of the dielectric constant, respectively, 
e’(w) 1s the frequency dependent value of the di- 
electric constant, ¢’’ the loss factor (€’’=e’ 
tan 5), tan 6 is the loss angle tangent or dissipation 
factor, f is the frequency in cycles per second, 
w=2nf, 7 is the relaxation time for dipole relaxation, 
n is the viscosity, @ is the radius of the spheres, k is 
the Boltzmann constant, and 7 the absolute tem- 
perature. 

Relation (2) has generally been found invalid, 
because in most systems investigated, values of the 
viscosity calculated from it differed from those 
observed, sometimes by orders of magnitude. The 
quantity » has therefore often been represented as 
an “internal viscosity,” a concept lacking quanti- 
tative definition in terms of measurable variables. 
In alternate treatments, applicable to systems with 
sharply specified energy states [4], the hydrodynamic 
theory is abandoned and a thermodynamic one is 
substituted. 

In the case at hand, it will be found useful to 
retain the hydrodynamic concept in a qualitative 
sense. This predicts a dependence of the relaxation 
time, and, therefore, of the frequency region in 
which dipolar losses occur, on viscosity and tem- 
perature. The usual way of studying the effect 
of these variables is, of course, by determining the 
effect of temperature on the dielectric properties 
of the solid polymer [5]. The hydrodynamic concept, 
however, suggests that similar effects might be 
observed by changing the viscosity at constant 
temperature. This picture, as we shall see, does 
indeed furnish a qualitative interpretation of the 
dielectric changes occurring during the polymeri- 
zation reaction. In fact, the agreement between 
corresponding changes occurring in the fully poly- 
merized polymer as a function of temperature and 
those occurring in the polymerizing mixture as a 
function of reaction time will be seen to be semi- 
quantitative. 


(w) 18 
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2. Materials and Methods 


2.1. Materials 


The primary ingredients of the casting resin 
were commercial grades of styrene, polystyrene, 
divinylbenzene, and hydrogenated terphenyl plasti- 
cizer. Of these the sytrene and the plasticizer were 
used without further purification, whereas the acrylo- 
nitrile (AN) was redistilled once, removing most, 
but not all, of the inhibitor. To insure a reasonable 
shelf life and a good curing rate of the stock solution, 
the composition of which will be described presently, 
the polystyrene used was a 50-50 mixture of poly- 
mers with weight-average-molecular weights of 
160,000 (manufacturer's data) and 840,000, as 
determined by light-scattering measurements. The 
divinylbenzene preparation (DVB) was stated by 


the manufacturer to contain 40 percent of divinyl- | 


benzene, most of which was the meta isomer. The 
other ingredients in this product were mainly diethyl- 
benzene and ethyl vinylbenzene. The inhibitor was 
removed from it by passing it through activated 
alumina. 

The minor ingredients of the resin were purified 
commercial benzoyl peroxide, a cobalt catalyst 
(cobalt naphthenate) containing 6 percent of cobalt, 
which was used without further purification and 
Promoter A.* 


2.2. Methods 


The systems studied, except the ones containing 
no DVB, were made up from a stock solution of the 
following composition by weight: Styrene, 50 parts; 
polystyrene, 30 parts; hydrogenated terphenyl, 13 
parts; AN, 5 parts; DVB solution, 2 parts. Systems 
termed 6-percent and 1l-percent DVB were made 
up by adding 4 g and 10 g of the DVB solution, 
respectively, to 100 g of the stock solution; the 
system with the nominal AN content of 10 percent 
was made up by adding an additional 5 g of AN to 
100 g of the stock solution. The system termed 
0-percent DVB was made up separately and differed 
from the stock solution in that the 2 g of DVB 
solution were replaced by an equal weight of styrene.’ 
Prior to each experiment, 0.4 percent of benzoyl 
peroxide, 0.04 percent of promoter A, and 0.15 


percent of cobalt catalyst were added to the other | 


ingredients to start the polymerization, and the 
viscous solution was poured into the cylindrical 


* This material is a 10-percent solution of B-hydroxy ethy! tetrahydroquino- 
line in hydrogenated terpheny! and is covered in a pending patent application 
by Max M. Lee 

* A casting resin made from a reaction mixture of this composition has been 
developed at N BS and termed A N-5 casting resin M. M. Leeand E. L. Hebb, 
unpublished experiments 

’ It was later found that aging of the stock solution may have had a slight effect 
on the experimental results. A fresh solution, and one which had stood at room 
temperature for several months, underwent electrical changes during polymeriza- 
tion in satisfactory agreement with one another, as far as the magnitudes of the 
electrical variables were concerned. However, the plots of these variables as a 
function of time were somewhat displaced toward shorter times in the case of the 
old solution, presumably because of a small amount of progressive poly merization 
This variation is of no importance, as far as most of the experimental results 
resented in section 3 are concerned, because we shall not be concerned with 

inetic questions, for the most part. In the one case where a qualitative examina- 
of kinetic factors was made (effect of peroxide concentration) care was taken to 
perform all the runs less than | month apart and to keep the stock solution under 
refrigeration. 








sample holder, which was kept at 50°+0.5° © 5 
circulating water from a thermostat through jj, 
jacket. The samples were kept at 50° C for aboyw 
35 hr and were subsequently cured for 14 hr at 109 
C. The 50° cure, during which the electrical seas. 
urements were carried out, took place in a room j 
which the humidity was controlled to not exceed 5 
percent, and the sample holder was kept stoppered 
during that period. 

Also, to make possible a determination of th 
dielectric properties of the fully polymerized resin as 
a function of temperature, the polymerizing mixtures 
were poured into forms somewhat larger than th 
required samples and were cured in a desiccator a 
the same temperatures and for the same lengths of 
time as the samples investigated during their 
polymerization. 

Measurements of dielectric constant and dielectr; 
loss were made throughout the polymerization 
two two-terminal sample holders designed for liquids 
and consisting of stainless-steel concentric cylinders 
separated by a \-in. gap. These sample holders 
were calibrated with benzene. Their interelectrod 
and residual capacitances were approximately 7 an 
2 uuf, respectively. Although this equipment should 
be inherently capable of dielectric-constant deter. 
minations accurate to several tenths of 1 percent 
some of the measurements probably resulted i 
values that were as much as 5 percent low, as 
discussed later. Loss tangents could be determine 
with an accuracy of +3X10~‘, resulting in los. 
factor accuracies of +9 10~‘ at a dielectric constant 
of 3.0. The same variables, as well as the d-c con- 
ductivity, were also measured as a function of tem. 
perature for the fully cured samples prepared in disk 
form. These measurements were made on a high- 
temperature two-terminal sample holder for solids 
housed in an air thermostat [6]. The estimate 
accuracy of the dielectric-constant determination o1 
this sample holder is believed to be within 1 percent 
that of the loss tangent was +2 10~‘*, resulting i 
loss-factor accuracies of +6X107~* at a dielectr 
constant of 3.0. The d-c conductivity measure- 
ments were made with the help of a d-c micr- 
microampere amplifier connected directly to th 
high-temperature sample holder, capable of measw- 
ing currents of 2X10-"\amp. The samples prepare 
for these measurements had been machined to giv: 
disks 1% in. in diameter and \, in. thick, will 
surfaces flat and parallel within 0.0005 in. Second: 
ary electrodes had been applied to these samples ) 
spraying them with a commercial silver paint 

All measurements of dielectric constant and los 
were carried out on the same model of a commereitl 
a-c bridge. Two different bridges, but of the sam 
make, were used with the cylindrical and the high 
temperature sample holders, which, because of 4 
replacement of a capacitor in one of them, gav' 
capacitance readings of different accuracies, i. ¢., 
+0.005 and +0.001 yuyf, respectively. The signal 
was generated by a beat frequency oscillator 
capable of a 10-v rms output. For detection of the 
signal a preamplifier and a cathode-ray oscilloscop 
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netion with the high-temperature sample holder. 
rogen analyses were performed by a micro- 
ahl method. These gave the acrylonitrile 
ntration in the cured polymer. Separate 
ses were performed for each sample investi- 
in the cylindrical sample holder, because small 
ences in acrylonitrile content were to be 
ted due to evaporation during handling and 

the polymerization. In the case of the 
nens prepared for the high-temperature sample 

only 1 sample with a nominal acrylonitrile 
nt of 5 percent and 1 of 10 percent were ana- 
because all these samples were cured in | 
and in the same desiccator, and because 
ences in acrylonitrile content between them 
be assumed to result only from differences in 
onitrile content in the polymerizing mixture, 
h can be calculated from the percentage com- 

position. 


3. Results and Discussion 


3.1. Electrical Behavior During the Polymerization 
Process 


Figure 1 illustrates the electrical behavior of a 
typical polymerizing mixture as a function of time at 
several frequencies. The initial factor, es- 
pecially at the lower frequencies, was very high and 
decreased rapidly with time. The frequency de- 
pendence of these initially high losses suggested the 
presence of considerable d-c conductivity, and this 
was borne out in a d-e conductivity measurement of a 
sample of acrylonitrile, which gave a value of 
approximately 2.5107? mho/em. This condue- 
tivity was much too high to be attributable to the 
acrylonitrile itself, and could be reduced roughly by 
a factor of 50 if a careful fractional distillation 
process Was carried out. This was not done in the 
experiments cited because the a-c losses attributable 
to the d-e conductivity decreased rapidly with time 
and vanished in the fully cured resin. The small 
amount of cobalt catalyst in the polymerizing mix- 
ture, because of its ionic nature, also contributed to 
he previously noted upswing of the loss curves 
toward the zero time axis. After these ionic losses 
had substantially disappeared, a second rise took 
n the loss factor at all frequencies, accom- 
panied by a decrease in the dielectric constant. This 
oss appeared and went through its maximum, first ai 
the highest frequency and last at the lowest frequency, 
ind is attributed to dipole orientation associated with 
arelaxation time, r. The loss factor must therefore 
co through a maximum at a frequency fax given by 
eq (1 In terms of the qualitative hydrodynamic 

the “internal viscosity”, and therefore the 
on time, will increase with the extent of the 
nm. This explains the observed shift of the 
oss Inaximum from high to low frequencies as the 
polymerization procedes. The gel point occurs 
shortly after the 70,000-c/s loss factor starts on its 
The loss factor shows no discon- 
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Fictre 1. Dielectric constant and loss factor at 560° Casa 
function of time after initiation of polymerization for a sam ple 
containing 4.9 percent of AN and 2 percent of DVB solution 
), 100¢/s; A, 1,000e/s . 70,000 c/s. Filled in symbols and dashed curves refer 


to an alternate run made with a polymerizing mixture of the same nominal com 
position for which an AN analysis was not available 


The curves of figure 1 bear 


tinuity at the gel point. 
a superficial resemblance to the Debye dispersion 
curves modified for a distribution of relaxation times 
(7, 8, 9] if the time axis is imagined replaced by a log 


frequency axis. This results from the fact that the 
rate of change of dielectric constant with time must 
be correlated with the magnitude of the out-of- 
phase component of the complex dielectric constant, 
i. e., the loss factor. Figure 1 also presents points 
from another run made with a fresher stock solution. 
The reproducibility is satisfactory as far as the shape 
of the curves is concerned; their displacement with 
respect to one another is probably caused by some 
progressive polymerization in the older solution, as 
discussed in section 2. A slight variation in the 
acrylonitrile content of the two samples may also 
have contributed to the shift. 

To a first approximation, the magnitudes of the 
loss-factor Maxima at any frequency are propor- 
tional to the concentration of acrylonitrile at a 
DVB concentration of 2 percent (table 1). The 
loss is therefore attributed primarily * to the orien- 
*E. B. Baker, R. P. Auty, and G. J. Ritenour, J. Chem. Phys. 21, 159 (1953 
have shown that polystyrene itself exhibits dipolar losses, which presumably 


result from the asymmetry introduced by the phenyl groups. This loss is 
inseparable from the one described, but amounts to only a small fraction of it 








TaABLe | Magnitude of loss-factor maxima during polymeri- | Taste 2. Time for 70,000-c/s loss factor to reach its m 
zation at 50° "C for various concentrations of AN and of DVB as a function of the peroxide concentration; polyme 
solution 7 carried out at 50° C 

Concentrations of AN and of DVB solution are 4.9 and 2 percent, res 

Concen 

tration a 
of DVB : . Peroxide 
) 70,000 c/s concen- rime 
tration 


solution 


wt. % 
0.1 


A similar dependence TaBLe 3. Effect of an additional 60-hour cure at 100° C 
the electrical properties at 24° C of a sample containing 


percent of AN and 2 percent of DVB solution 


tation of the CN dipoles. 
the magnitudes of the loss-factor maxima on the 
acrylonitrile concentration is found for the high- 
temperature measurements made on the fully cured 


Standard cure Prolonged cure 


polymer discussed later (figs. 4 and 7, table 5). 

The peak losses at any frequency decrease as the | enite 
concentration of cross-linking agent is increased , eee oe 
(table 1). A plot of loss factor versus log frequency on 
for various concentrations of DVB shows the maxima | c/s ; : 
to decrease and the curves to become broader as 500 2. 5 O08 ‘ 0. 0020 
the concentration of DVB is increased (fig. 2).° bei 3 — - = 
These phenomena can best be described in terms of | 70. 000 25 og | 25 * 0020 
a broadened distribution of relaxation times result- 
ing from an increasing failure of the CN dipoles to 
orient independently of one another. This effect is 







































seen to be most pronounced at the highest DVB | from 0.2 to 0.4 percent. This indicates that the tiny hie} 

concentrations. | at which the loss peak occurs at any one frequen — 

Table 2 shows the time at which they 70,000-c/s | under conditions where the structure of the polym er. 

loss peak occurs for various peroxide concentrations. | and the reaction mechanism may be assumed to } “ " 

Qualitative observations of the fluidity of the solu- | identical, is a convenient measure of the reaction rat “a ‘ 
tion and, later, of the rigidity of the gel as a function After approximately 35 hr at 50° C, followed by 4 

of time indicated that the polymerization rate was | 14-hr cure at 100°C, this being the standard ew erp 

strongly dependent on the peroxide concentration at | condition, the dielectric constant and loss at all fr- — 

concentrations below 0.2 percent, but that no such | quencies had become essentially constant. Table in a 

dependence could be readily observed in the range | shows the effect of an additional 60-hr cure at100° | . Ms 

There is no significant change in the values of the los a 

factor during that period. The new values of tl a 

rs T dielectric constant are slightly lower, but this is pro)- a 

a a ably caused by shrinkage of the sample. The concl- i , fo 

a ig sion seems warranted that the electrical properties ce ni 

5 .2s- the sample did not change significantly after th atta; 

Sm ' initial cure. This conclusion is not surprising, as, elects 

ak terms of the qualitative picture presented here, a no- ar 

polar plasticizer and small amounts of unpolymerize ne - 

monomer would be expected to affect the electric Surf, 

properties in a similar manner, i. e., primarily mer 

a reducing the “internal viscosity.”’ It has been show! <a 

_ by Fuoss [10] that the effect of the plasticizer is mos with | 


* 9 7 , . ; 5 ’ > > th 
vom - : noes factor — oe her at my p. = pronounced at very low concentrations. In th for tl 
Junction of frequency a e approrimate ime at which the | — . a = : stimiwar 7 _ . ys ¢ . 
1 000—c/s loss-factor marimum occurs for several concentrations poor nce of “ latively high plastic ser comsenires 
of the DVB solution (13 percent in this case) small amounts of unpoly- 
Loss-tactor values have been corrected to an AN content of 5.0 percent merized monomer should therefore have no signi! 
xs8-facto s hav » ~~ “ter >a content of 5 vercen . . ~ } f 
percent DVB solution. A, 2-percent DVB solution Coureunt pvB. | cant effect on the electrical properties, which, | 
. tio yr, il Tce VB so o . : } 
on pact DVS eteien | systems like the one discussed here, should not! 
expected to provide a sensitive criterion for col- 
on In qrder to make this comparison more significant, the data plotted in this pleteness of cure. 
gure have been corrected to the same AN concentration. This is done by rv. ’ i ° artes 
making use of the fact that any linear theory of polarization requires A. Table 4 prov ides data A omparing the prope - crepal 
tionality between dipole concentration and dielectric loss at any frequency, | Of a sample cured as described and measured in th 


subject to the restriction that the distribution of relaxation times is independent . . ° : . t one ni 
of dipole concentration cylindrical cell for liquids with the average comput twe b 
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Comparison of electrical properties of two samples of 
I composition, c ured in cylindric wu cell and disk forn . 
vely 


Sample A 


Dielec- 
ency tric con- 
Stant 


Loss 


factor 


Loss 


er tri 
factor Frequen ri 


y 


c/s 
0. 0026 100 
0030 1, 000 
0028 10, 000 
OO1L6 100, 000 


0. 105 loo 0. OS9 
O78 1, 000 4 O55 
032 10, 000 > & 033 
O19 100, 000 > 025 


e cylindrical cell. (Same sample used in table 3 
mputed from ? samples cast and cured in disk form and measured 
temperature holder. The AN concentration of this sample was 
of sample A 


two fully-cured disks of identical composition 
complete data are presented in fig. 4) measured on a 
high-temperature sample holder. Data for this 
comparison had been made available by taking 
several of the samples polymerized in the cylinderical 
cell through part of the temperature range covered 
by the high-temperature sample holder and by in- 
terpolating (fig. 4) from values obtained on that 
sample holder to make the comparison at 80° C. 
In accounting for discrepancies between the two sets 
of values, the following factors should be kept in 
mind: Shrinkage of the sample during polymeriza- 
tion in the cylinderical cell, curing conditions that 
could not be made exactly identical (e. g., different 
surface areas and different rates of evaporation), 
different geometries and different electrical arrange- 
ments of the sample holders, and different rates of 
attainment of temperature equilibrium throughout 
electrodes and sample. Shrinkage during polymeri- 
zation probably accounts for the lower dielectric- 
constant values obtained in the cylinderical cell. 
Surface conductance at the lower frequencies, and a 
small effective series resistance due to the secondary 
silver electrodes at the higher frequencies, together 
with the difference in curing conditions could account 
for the somewhat higher loss-factor values obtained 
on the high-temperature sample holder at room 
temperature. The discrepancies between the loss- 
factor values obtained at 80° C are probably due to 
imperfect temperature equilibrium, because in this 
region the loss factor changes rapidly with tempera- 
ture. The most rapid change takes place at pre- 
cisely the frequency (1,000 c/s) at which the dis- 
crepancy is largest. In view of these considerations, 
one may conclude that the samples prepared by the 
two methods were closely similar in their electrical 
Tespor It should be kept in mind that studies of 


dielectric properties as a function of temperature are 
as yet relatively rare and that questions of the 
reproducibility of such measurements at elevated 
temperatures, even when carried out on the same 
sample holder, have been left virtually unexplored. 


3.2. Electrical Behavior of the Fully Cured 


Resin 


When the electrical properties of the fully cured 
resin are investigated as a function of temperature, 
a familiar behavior is found [5]. At any frequency 
there is a sharp rise in dielectric constant with 
temperature, a maximum is reached, followed by a 
slow decrease. Paralleling this, there is a rise in loss 
factor, a maximum occurs, and by the time the 
dielectric constant has gone through a maximum, 
the loss factor is small again. It is also a familiar 
fact that the rise in dielectric constant and loss 
factor shifts to higher temperatures as the frequency 
is raised. This is consistent with eq (1) and the 
qualitative hydrodynamic concept, which predicts 
a decrease of the relaxation time with temperature, 
keeping in mind that + represents an average value 
introduced because we are dealing with a distribution 
of relaxation times. The behavior described is 
illustrated in figures 3, 4, 5, 6, and 7. Figure 4 
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| Fiaure 3. Dielectric constant and loss factor of the fully cured 
resin as a function of temperature for a sample containing 
O percent of the DVB solution. 


AN content is 5.6 percent , 00 c/s, A, 1,000 c/s , 10,000 c/s; WV, 100,000 c/s, 
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Fiaure 4. Dielectric constant and loss factor of the fully cured 


resin as a function of temperature for a sample containing 
Zz pe recent of the DVB solution 


AN content is 5.6 percent 100 e's 1,000 cs 1,000 ¢ 100,000 c/s, 


Filled in symbols refer to an alternate sample 


includes points from a duplicate run. As the 2 
samples from which the 2 runs were made were 
prepared under identical conditions, the discrepancies 
between them are probably largely attributable to 
imperfect temperature equilibrium. The initial small 
decrease in dielectric constant at the higher frequen- 
cies would appear to indicate that in this region the 
disordering effect of temperature overshadows the 
dipole alinement. 

Figure 8 shows the loss factors as a function of 
frequency for several concentrations of DVB at 
approximately the temperature at which the 1,000-c/s 
loss factors go through their maxima. We note here 
a parallel with the behavior of the resins during their 
weg rman at the times at which the 1,000-c/s 
oss factors go through their maxima in those systems 
(fig. 2). In both cases we find that the 
generally become less and the loss factor versus log 
frequency curves become broader as the concentra- 
tion of DVB is increased. The interpretation in 
both cases is the same, although for the fully cured 
= the effect only becomes definitely estab- 
ished at the highest concentrations of DVB. We 
also note that the magnitudes of the losses in both 


losses 





DIELECTRIC CONSTANT, € 








LOSS FACTOR, ¢” 





L i i i i 
20 40 60 80 100 120 140 
TEMPERATURE, °C 





Ficure 5. Dielectric constant and loss factor of the fully 
resin as a function of te mperature for a sample containing ' 
percent of the DVB solution. 

, OOes: A, 1000 e's , 10,000 e's 


AN content is 5.3 percent. 


cases are similar. This can be seen best by con- 
paring the loss maxima for each DVB concentration 
and for each frequency, that is, by picking the proper 
time for the polymerization studies and the proper 
temperature for the high-temperature studies on thi 
fully cured polymers (table 5). This observatior 
can be rationalized by the supposition that |) 
losses depend essentially on the average dipol 
relaxation time, that changes in this relaxation tim 
can be brought about by variation of temperature or 
by changes in the “internal viscosity’? occurring 
during the polymerization at constant temperatur 
and that the distribution of relaxation times is no! 
greatly different in the systems compared. Becaus 
average molecular weights and molecular-weigh 
distributions during and after complete polymeriz- 
tion are different, this would also mean that th 
relaxation distribution function is essentially inde 
pendent of these variables throughout the range cot- 
sidered. It is worth noting that these systems 0 
equivalent electrical response are rather different 
their macroscopic mechanical appearance—the poly- 
merizing systems can be described best as gels, th 
fully polymerized systems as rubbers. 

The temperature of the loss-factor maximum # 
any frequency increases with the concentration © 
DVB. This is related to the observation that the 
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LOSS FACTOR , €" 





Ficure 8. Loss Jactor of the fully cured resin as a function of 
frequency at the approximate temperature at which the 1,000- 
c/s loss-factor marimum 


the DVB solution 


several concentrations of 


occurs for 


AN content of 


Loss factor values are corrected to an 
olution 


DVB solution 2-percent DVB 
1l-percent DVB solution. 


».0 percent 0-percent 
t-percent DVB solution 


TABLE 5 Comparison of loss-factor maxima reduced to AN 
contents of 5 or 10 percent of polymerizing systems at 50° ( 
and of fully polymerized systems at elevated tem peratures 


100,000 « 


5, 6, and 7 


loss maxima occur in the temperature range in which 
the samples undergo the change from glass-like to 
rubber-like behavior [11, 12], and that this transition 
occurs at higher temperatures the higher the con- 
centration of DVB. 

At the higher temperatures the samples exhibit 


some d-c conductivity (table 6), which could be due 
to a small amount of pyrolysis or to a freeing of 
ionic impurities in the polymerizing mixture. lonic 
conductance in phase with the applied field gives 
rise to an equivalent a-c loss factor equal to 1810" 
G/f, where G is the d- conductivity. When this 
quantity is subtracted from the measured a-c loss 
factor, the remainder should be representative of the 
dipolar loss alone. A calculation shows that after 
this subtraction is made, the a-c loss factor still 
increases with temperature, after having gone 
through a minimum, and, although this may mean 
that a second dipolar loss mechanism has come into 
play, an alternate explanation is that the ionic 





TABLE 6 Direct-cu 


0 percent of DVB solu 6 percent of DVB solu- 
tion, 5.3 percent of AN 


{ DVB solu 
5.6 percent of AN reent 


tion {AN 
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ictivity 
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ent conductivities of fully cured resins as a function of te m perature 
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10-44 


conductance has an out-of-phase component, giving 
a larger a-c loss than that which would be expected 
from a pure ionic in-phase mechanism [13]. 

By application of the theory of absolute reaction 
rates to dielectric relaxation [14], one obtains 


kT /he®5/*®¥e-48/8T (3) 
where & is the Boltzmann constant, A is Planck’s con- 
stant, and AS and AH/ are entropy and heat of activa- 
tion, respectively. Calculations of AS and AH were 
performed for all samples for which sufficient data were 
available (all except the sample containing no cross- 
linking agent). As an additional check, a run was 
made with another sample containing 2 percent of 
the DVB solution, and measurements were made at 
closer frequency intervals (table 7). Figure 9 shows 
the plots of log fas. versus 1/7’ for all these systems 
made from the data in figures 4, 5, 6, 7 and table 7. 
These plots give, within the experimental error, 
parallel lines, and the heat of activation averages 
63,000 cal/mole. For the entropies of activation, 
values of 122, 115, and 107 (cal/mole)/deg are ob- 
tained for the samples containing 2, 6, and 11 percent 
of the DVB solution, respectively, and the normal 
amount of acrylonitrile, whereas a value of 115 
(cal/mole)/deg is obtained for the sample with the 
high AN content. The systems investigated there- 
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| TABLE 7. Detailed high-te mperature data for a samp: 
taining 9.4 percent of AN and 2 percent of DVB sea 


Loss 
factor 


Fre- 
quency 


Dielectric 
constant 


Loss 
factor 


Dielectric 
constant 


Fre- 
quency 


108° C 





fore show the high values of both heat and entro; 
of activation commonly observed for dipolar relax- 
tion in high polymers [11]. 

It has been shown that distributions of relaxatio: 
times can be calculated from dielectric dispersio 
[9, 15, 16] and that dielectric losses can be calculat 
from dielectric-constant data without introducing 
any assumptions concerning the form of the distn- 
bution function [15]. We have for a linear dieclectr 
of very small d-c conductivity: 

= y(r)dr 


, 
€ (w)—e€, i >= 
J 0 1+ w*r* = 
. . 
“- “y(r)wrdr . 
€ (a) = - —— = 
0 l tT @ Tt” on 


ao 


+? Yiln r)dlnz 
= 1+w’*r’ 


wT Y(n r)d \n-: 


1+ 3*r’ 


and, if the dispersion is not too sharp, we have to: 
very good approximation [15): 


de’ 


Y (In 1/e)= 7 





where y(r) is the distribution of electrical relaxatios 
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and Y(In r)d In r= y(r)d7 Often it is of 
to obtain the function Y from*the imaginary 
well. This has been attempted for the analo- 
echanical case [17, 18]. The same treatment 
| to the electrical case would result in 
op d log ” ad w) 


Y (In 1/w) € 


d log Ww 


been pointed out [15], the use of the approxi- 
n resulting in (7), when used in the calculation 
rom e’, leads to values that are 30 to 40 percent 
iall. The reason for this discrepancy is that, 
is the approximation made in the calculation 
namely, 


(1 0- 
o 


ellent, the approximation which results in 


ery 
{or UV. 


lo 


because we are interested in the integrals of 
expressions (fig. 10). It is clear, therefore, that 
cannot obtain Y with the aid of approximation 
We can, however, define a function Y’, given by 


* In l/w 
ee | wrY’(Inr)dinr- (10) 


— 2 


not r. 
which 


this funetion Y’, 
hands side of (7) 
aced by 


that is equal to the 
must, therefore, be 


nitro} 


relaxa- ” d log (e’’ wW) 


Y’(in 1 w) *g (11) 


d log WwW 
xatiol 
eTSIO! 
ulate 
ducing 
distr- wrY(Inr)din + 


‘leet 1+w*r* 


can now calculate Y’ 
0), obtaining 


in terms of Y from eq (5) 
ein 1/w 
wrY’(Inr)dIin-r- (12) 


wines 


and y’ 


in most practical cases, the variation of Y 
is small in the region in which the integral 
ignificant value, we can write instead of (12) 


48 f * inl /w 
Y (In 7 | 7 dinr= “dn r) | wrdinr- 
» i+w 


re) 


(13) 


Making use of the symmetry of the integral on the 
leit, and transforming, we obtain 
*l/e *l/e 
9 WwW , 
2) (In 7) —, dr=Y'(Inr) wir- (14) 
0 l TW T* 0 


uxation Pr 


speaking, this is a distribution of relative permittivities. 




















-2 
WT 


Fiaure 10 and of 
required for correlations between real and imaginary parts of 


») . = 
Plot of exact €L Pressions approximations 


the dielectric constant plotted as a function of log wr 


On the left, 1/(1+*r? 
approximation (9 


and approximation (8 
exact expression 


on the right, wr! (1+ and 
approximations, 


Hence 
Y’(In + 


r/2 Vile r (15) 
i. e., the two functions are simply related through a 
constant over their entire spectrum, and by (6), 
(11), and (15) 


de’ 
—1/2.303 —? 
d log f oil 


P d log Pf / 


Y (In r) 
d log J 


16) 


The appearance of the factor 2/r immediately 
explains the discrepancy between the observed and 
calculated values of ¢’’ noted previously [15], since 
the use of approximation (9) is equivalent to the 
neglect of the difference between Y and Y’ For 
materials whose losses are essentially frequency 
independent (12) reduces to 


. de’ 
Y(In r)=—1/2.303 , 
d log / 

The identical relation has been derived by Gevers 
and Du Pre, using a somewhat different argument [19]. 
Equation (17) provides a very convenient first 
approximation connecting dielectric dispersion and 
Table 8 gives a comparison of values of the 
loss factor as observed and as calculated by the 
procedure previously described [15], i. e., through 
use of eq (6), followed by (5). It also presents 
values of )(In +) calculated from real and imaginary 
parts of the dielectric constant through the use of 
eq (16). The reasonable agreement of the two 
sets of e«’’ and }) values demonstrates the self- 
consistency of the approximations outlined, which 
is fully adequate in view of the inaccuracies en- 
tailed by the graphical operations. 

For the Saalanees mechanical case, successful 
calculations of dynamical properties at various tem- 
peratures have been made from measurements at a 


Cm 


kc SS. 


"Ferry et al. [17], who use approximation (9), noted that in some of their 
results the use of an empirical constant ‘‘a’’ equal to 0.77 occurring in the place 
of 2/r in an equation analogous to (16) harmonized distribution functions calcu- 
lated from real and imaginary parts of the modulus 
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TABLE 8 
percent of AN and 2 percent of DVB solution at 94° ¢ 


ob e” (cal Y (from | Y (from 


Fre e : 
culated 


quency served « « 


P22 PS tO Pe PO PO PO te we 


different temperature, using the assumption that the 
distribution of relaxation times did not change with 
temperature [20]. Although some success has been 
reported for the electrical case also [16], similar cal- 
culations of ¢’ and ¢’’ as a function of temperature, 
from data obtained at a different temperature, do 
not give sufficiently accurate results when applied to 
the data cited here. The reason seems to be that 
the distribution functions are not independent of 
temperature. As it is not possible to calculate dis- 
tribution functions with an accuracy of better than 
5 to 10 percent by the method outlined, we may cite 
the loss-factor maxima at 
as evidence that 
Since 


the observed increase in 


the higher temperatures (table 5 
these functions sharpen with temperature 


@—€,. iS proportional to e’’d In f, and since & 


. 
decreases with temperature while ¢, is essentially 
temperature independent over a small temperature 
range, an increase in the loss-factor maximum with 
temperature must be connected with a sharpening of 
the distribution function with temperature.” Such 
a sharpening has been concluded to occur for 
polymers in general by several investigators, both on 
the basis of experimental data and on theoretical 
grounds [9, 11, 12] 

Because of the analogous effects on the relaxation 
time of temperature changes in the case of the cured 
polymer and of “internal viscosity’? changes in the 
case of the growing polymer, it is interesting to note 
that the loss-factor maxima during the polymeriza- 
tion increase with frequency (table | 


4. Summary 


1. The marked changes in dielectric constant and 
in dielectric loss observed during the polymerization 
of the polar casting resins described are attributable 
primarily to the dipolar relaxation of the nitrile 
groups present 

2. As seen by its electrical effects, the hardening 
of the reaction mixture at a certain stage of the 
polymerization is quite similar to the glass-phase 
transition occurring when the temperature of the 
fully cured rubber-like polymer is lowered—the 


4” This argument assumes that the resultant dipole moment of the polymer 


molecule is independent of temperature. 


Electrical properties of a sample containing 5.4 | 


dipolar relaxations occurring in the above typ 
transitions result in loss factors that agree sem). 
quantitatively, and a similar broadening of ¢h, 
relaxation distribution functions take place ipo; 
addition of cross-links in each case. 

3. In the polymerizing systems discussed, ; 
trical measurements can be used as an indicator oj 
reaction rate in a qualitative way. However 
questions as completeness of cure cannot be decided 
in this manner because all significant electrica 
changes occur at a relatively early stage of the 
reaction. 

4. A scheme is described whereby distributions of 
electrical relaxation times can be calculated wit} 
comparable accuracy from real and imaginary parts 
of the complex dielectric constant. , 


suc 
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Study of Degradation of Polystyrene, Using 
Ultraviolet Spectrophotometry 


Mary Jane Reiney, Max Tryon, and B. G. Achhammer 


The role of ultraviolet radiant energy and effect. of presence of monomer on degradation 


of polystyrene were studied by means of ultaviolet absorption techniques 


It was found that 


exposure of pure polystyrene to ultraviolet radiant energy resulted in increased absorption 


in the ultraviolet reg 
time increased 
ing the general tv pe ot absorption 


polymer and polymer containing monomer styrene 
crease d the rate of progress of the post-radiation effect 


on of 280 to 400 millimicrons, which advanced progressively as exposure 
rhe presence of monomer increased the rate of degradation without alter- 
\ post-radiation effect was noted in both the purified 


Again, the presence of monomer in 
Possible mechanisms are postulated 


for the ultraviolet degradation, and the concept of entrapped free radicals is considered as a 


possible explanatior 


l. Introduction 


s been shown by infrared spectrophotometric 
ments that carbonyl and hydroxyl groups 
med in polystyrene on exposure to ultraviolet 
energy in the presence of oxygen [1, 2].' 
latile products given off and measured by 
ectrometry and the changes in solubility and 
content on exposure of polysty rene to heat 
iltraviolet radiant energy in vacuum and in 
have also been reported [3]. This report 
the changes in the ultraviolet spectrum of 
ene as a result of exposure to ultraviolet 
energy in air. A mechanism is suggested 
ribe possible chemical changes in the polymer 
ould be responsible for the variation of the 
let spectrum of polysty rene, both on exposure 
polymer to ultraviolet radiant energy and on 
ent storage in the absence of light. 


2. Materials and Methods 


‘he polystyrene samples used in this study are 
scribed in table 1. All samples were prepared by 
Dow Chemical Co. Samples A and El are 
mercial resins; they do not contain any com- 
unding ingredients and were used as received. 
samples W, X, and Y were prepared by polymerizing 
irified monomer without catalyst or solvent in a 
trogen atmosphere. 


ABLE 1 Descr iption oj sam ple 8 of polyst yrene 


Approximate Polymeriza- 
molecular tion tem 
weight perature 


Sample 


A * 308, 000 
El * 200, 000 
W > 500, 000 
x b 237, 000 
Y > 110, 000 


* Weight average obtained by light scattering. 
> Number average obtained from osmotic pressure 


— 


rackets indicate the literature references at the end of this paper. 


for the post-radiation effect observed 


The spectrophotometric measurements on these 
samples as received are designated by the sample 
letters alone. The measurements made after partial 
purification of the samples by a “frozen benzene 
technique” [2] are designated by the sample letters 
followed by P. Sample X was also subjected to 
purification by solution in methyl ethyl ketone and 
precipitation with isopropanol three times; the 
samples from the first, second, and third precipita- 
tions are designated X1, X2, and X38, respectively. 

The technique for preparation of the samples in 
film form of various thicknesses has been described 
previously |2]. Casting of the thicker films on a 
surface of regenerated cellulose, instead of directly 
on glass, facilitated their removal because wetting 
the opposite side of the cellulose with water swelled 
the cellulose and broke the bond between it and the 
polystyrene film. 

Exposure of the films to ultraviolet radiant energy 
was accomplished by placing the films in holders on a 
revolving table 6 in. from an S-1 sunlamp. The 
exposure was done in air, and the temperature at 
the table was 60° C. This sunlamp equipment is 
described in method No. 6021 of Federal Specification 
L—P—406a [4]. 

A Beckman model DU spectrophotometer with 
ultraviolet accessories was used for all the transmis- 
sion measurements.’ For solution work, silica cells 
of 1.000-cm path length were used. 


3. Effect of Film Thickness on Transmittance 


In the same manner that the concentration of a 
solution for a given light path length will dictate the 
spectrum of polystyrene obtained, the thickness of 
films of the polymer will be the controlling factor in 
determining its absorption spectrum. The effect of 
film thickness on the spectrum of polystyrene films 
of sample XP is shown in figure 1. 

Films of thickness equivalent to that of curve E 
of figure 1 were selected for the present study because 
they are thick enough to be mechanically stable; 

? The various spectral transmission and absorption terms used in this paper 


are based on the terminology recommended by Kasson 8. Gibson in N BS Circular 
484, Spectrophotometry (200 to 1,000 millimicrons) (Sept. 15, 1949 
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~ FULMS OF PARTIALLY PURIFIED SAMPLE xP) 


a in 0.00004 IwCH 


— 


28 ux 
WAVELENGTH, MILLIMICRONS 


Effect of film thicknesa o ultraviolet 


tlance of pr / 


transmi 


they are not too thick to show the shape of the 
absorption curve in the region 275 to 350 my where 
the differences between the polymers 
studied were obtained : and their strong absorptance 

below 270 mu is relatively unimportant to this study 

Smakula [5] states that the absorptance of poly- 
stvrene at about 260 ny 1s due to the phenyl residues 
These would not be expected to change appreciably 
on degradation of the polymer. Furthermore, mono- 
mer styrene, which has been suggested as playing an 
important part in the discoloration phase of the 
degradation of the polymer [1], is known to absorb 
in the ultraviolet at 291.5 mu. The fact that the 
absorption band at 291.5 mu is present in curve E 
indicates that the purification of the sample was not 
complete and that some monomer remained in the 


greatest 


sample 


4. Transmittance of Undegraded Polymers 


Figure 2 presents the spectra of films of the various 
polymer samples studied. It is thought that the 
differing spectra cannot be ascribed to the moderate 
thickness variations of the films. This is substan- 
tiated by the measurements made on solutions of 18 
g/liter concentration of WP, XP, and YP in benzene, 
which are presented in figure 3. The solution spectra 
of these three samples bear out the relationship shown 
by the spectra of these samples in film form. 

Moreover, the variation in the spectra cannot be 
ascribed to the molecular weight differences of the 
samples. McGovern and coworkers [6] found that 


the same spectrum was obtained at wavelengths of | 


282 and 291 mu, regardless of the molecular weight 
of the polystyrene. 
specific extinction coefficient, referred to weight con- 
centration of polystyrene, at 260 my is independent 
of molecular weight. Newall [8] measured five poly- 
styrene samples and observed that neither molecular 
weight nor variables of the polymerization reaction 


+ Absorptance is defined as “1— Transmittance." 


Meehan [7] showed that the | 


TRANSMIT TANCE 


SAMPLE MOLECULAR WEIGHT  FiLa 
ve 200 
t 293,000 
a 308,000 
mp 237,000 
we $00,000 


i 
320 34 
WAVELENGTH, MILLIMICRONS 
Ultraviolet transmittance 
purified pol ystyrenes 


of partial 
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SOLUTIONS ARE i6g/] IN GENZENE 
SAMPLE MOLECULAR WEIGHT 
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~ 237,000 
we $00,000 
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Ficure 3. Ultraviolet transmittance of partiall 


purified porystyrenes 


affected the spectrum of polystyrene at the shorter 
wavelengths. 

The spectrum of polystyrene, however, is a fune- 
tion of the purity of the sample. Sample X was care- 
fully purified by repeated solution in methy! ethy! 
ketone and precipitation with isopropanol. The \ 
curve in figure 4 represents the pure polymer, avd 
the other spectra denote varying degrees of contami 
nation retained by the polymer. The X polymer show: 
marked regions of absorptance not characteristic 0! 
| pure polystyrene. Figure 5 combines the spectra 0! 

18 giliter solutions of X3 and X with those of th 

unpurified polymers W and Y. The difference 
| purity of the polymers W, X, and Y is clearly show 
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Effect of monomer and impurities on the ultraviolet 


transmitlance of polystyrene. 


larly the considerable contamination in the 
W polymer. It is reasonable, therefore, that an 
dentical purification procedure applied to each poly- 
mer might not produce samples of the same degree of 
purity, thus leading to the differing spectra exhibited 
by the samples WP, XP, and YP in figure 3. 


5. Effect of Styrene Monomer on Polymer 
Spectra 


ln order to determine the part that residual mono- 
mer contributes to the contamination of the unpuri- 
led polymer, styrene was added to solutions of pure 
polystyrene (precipitate X3 of sample X). Solutions 
containing approximately 0.1, 0.4, 0.8, 1.17, 1.6, 
2.0, 2.5, 5.0, and 7.35 percent of styrene were pre- 
pared. Figure 6 shows the curve of the pure polymer, 
those of the first five solutions in the monomer- 


TIONS ARE 169// IN CHLOROFORM 
THIRD PRECIPITATE OF x 
INPURIFIED 180°C SAMPLE 
UNPURIFIED 20°C SAMPLE 
UNPURIFIED 60°C SAMPLE 
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ultraviolet transmittance of pol rslurene 


addition series and that of the unpurified polymer, 
all at 18 g/liter concentration. The effect of in- 
creasing monomer content is clearly shown. There is 
also an absorption band from 310 to 385 my in the 
unpurified sample. That this band cannot be 
attributed to monomer contamination is evident 
from the spectra in figure 7 made with solutions ot 
9 g/liter concentration. Even when the 292 my 
monomer band of X is approximately duplicated by 
the addition of 2.8 percent of monomer to X3 there 
is still no indication of the 310 to 385 mu band of X 
present in the spectrum of the latter solution. 
Changes in absorptance in the 310- to 385-myu wave- 
length region were neted on exposure of the films to 
ultraviolet radiant energy, and a more detailed 
investigation of this region is described in section 8. 

This failure of residual monomer to account for 
the long-wavelength contamination of the unpurified 
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polymers is further shown by examination of figure 
8. The pure polymer solution containing 7.59 per- 
cent of added monomer shows a spectrum com- 
parable to that of W only below approximately 
205 mu. The presence of monomer causes ho ap- 
preciable change in the spectrum of the pure poly- 
mer at wavelengths longer than 300 my. This figure 
shows the two monomer absorption bands at approxi- 
mately 282 and 291.5 mu more clearly than any of the 
preceding figures and illustrates that monomer 1s 
responsible for the decreased transmittance * of 
W below 295 mu 

Figure 9 presents the spectra of three solutions of 
stvrene monomer alone in chloroform. The absorp- 
tion bands at 282 and 291.5 my are shown, as well 
as the intense band at 248 mu; but even when the 
spectrum of a very concentrated solution is con- 
sidered (curve C), there is no indication of an 
absorption band at longer wavelengths. The absorp- 
tion bands shown are due to the presence of the 
double bond in conjugation with the benzene ring 
and form the very characteristic spectrum that ts 
distinguishable from that of saturated derivatives 
and from that of unsaturated compounds without the 
double bond in conjugation with the benzene ring [9] 
This pattern of absorptance displayed by the mon- 
omer permits quantitative determination of stvrene 
in poly styrene spectrophotometri ally 6, S| as well 
as the determination of monomer styrene in ethyl- 
benzene [10] 


6. Effect of Monomer on Polystyrene 
Degradation 


To study the effect of styrene monomer on the 
degradation of polystyrene, two films were pre- 
pared. One was cast from a benzene solution of 
sample X3 and the other from a benzene solution 
of sample X38 to which 2 percent of monomer was 
added. The films were made by the same procedure 
and measured 0.0066 and 0.0068 in. thick, respec- 
tively. They were exposed to ultraviolet radiant 
energy from an S-1 sunlamp at 60° Cin air. Ultra- 
violet transmission measurements of these films at 
intervals during their exposure give an insight into 
the degradation of the pure polymer and polymer 
contaminated with monomer but not with other 
impurities. Figure 10 shows the effect of exposure 
to S-1 radiant energy on the transmittance of the 
pure film in the ultraviolet region of the spectrum, 
and figure 11 shows similar information for the film 
prepared from the solution to which monomer was 
added. A progressive general increase in absorp- 
tance is noted in each case but is initially more 
severe for the film containing monomer. The 20- 
and 50-hr spectra of figure 11 bear a resemblance to 
the partially purified XP film spectrum noted 
previously. Both films yellowed very slightly as 
they were exposed, the monomer-containing film 
slightly more so than the pure film 


‘“Transmittance,”’ 7’, is defined as the ratio ot the radiant energy leaving the 
sample to the incident radiant cnergy; 1007 =percentage transmittance 
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7. Post-radiation Changes in Folysiyrene 
Degradation 


In the course of the degradation studies. a post- 
radiation effect was observed. The first indication 
of this effect was noted when, following the 15-hr 
S-1 sunlamp exposure, the films were stored in the 
dark at room conditions in air for 99 br and their 
transmittance measured before further exposure 
The resulting spectra showed a noticeable increase 
in absorptance when compared to the 
measured directly after exposure to radiant energy 
The 20- and 25-hr total S | exposures followed the 
pattern of the 15-hr exposure. After 155 hr of total 
S-1 sunlamp exposure, the films were stored for 13 
days, and a very large change in absorptance was 
apparent when the spectra were measured. The 
films were then exposed to an additional 5 hr of 
radiant energy from an S—-1 sunlamp, to give a total 
of 160 hr of exposure, and their spectra remeasured 
The increased absorptance apparent after storage 
was eliminated, and the 160-hr spectra fell in line 
with the spectra measured immediately upon com- 
pletion of the 155 hr of exposure. This effect was 
checked by storing the films for 13-day periods 
following the 200-, 265-, and 290-hr total exposure 
times and was found to be real. Figures 12 and 13 
present the spectra of the films measured just before 
and after storage periods. The dashed lines repre- 
sent the measurements made following the storage 
periods 

Upon completion of the 290-hr S-1 sunlamp 
exposure, the spectra of the films were measured at 
intervals during the ensuing storage of the films in 
order to determine the length of time necessary for 
the development of the increased absorptance. When 
the corrected absorbance, A,,° at a particular wave- 
length is plotted versus the hours of storage following 
the 290-hr exposure to ultraviolet radiant energy, a 
curve is obtained for each film which shows that 
most of the increased absorbance develops in the 
first 75 hr of storage. Figure 14 is such a plot for 
numerous wavelengths. From this figure it is seen 
that the slopes of the curves are very slight at 
greater than 75 hr of storage and also that a greater 
and more rapid change occurs at any particular 
wavelength on storing the monomer-containing film 
than on storing the pure film 

From this information, it is clear that the first 
storage-period measurement, which was made after 
99 hr of storage following the 15 hr of exposure to 
ultraviolet radiant energy, is representative of the 
maximum increase of absorbance due to storage 
following that S-1 exposure. Two hours of exposure 
to ultraviolet radiant energy is sufficient to eliminate 
the absorption increase due to storage. The 202-hr 
curves and the 267-hr curves were not presented in 
figures 10 and 11 because they closely follow the 
200- and 265-hr curves, with slightly decreased 
transmittance One hour of ultraviolet radiant 


—————_—$ ————-—- 


spectra 


§ 4. represents the absorbance reading of the film at a given wavelength minus 
the absorbance reading at that wavelength measured on the film prior to all 
treatment. “Absorbance”’ is defined as logy Js/J; also known as optical density 
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energy, however, was shown to be inadequate to 
cause this effect, for the 266-hr total-exposure curves 
fell below 267-hr curves and were of a shape more 
like that of the storage curves 

Because 2 hr of exposure to ultraviolet radiant 
energy at 60° C in air would return a storage curve 
to the previous 5-1 exposure curve, 2 hr of exposure 
60° C heat in air was administered to the films 
following the storage after the 290-hr total S-—1 
exposure to observe the effect of heat alone. Heat 
did not remove the storage absorption band, and 


there was no change in the transmittance as a result 
of the treatment 

\ graphical presentation of which 
depicts the complete pattern of treatment of the 
films, is shown in figure | (As representative of 
the absorption band covering the 310- to 385-my 
wavelength region, transmittance at 340 my 
plotted versus exposure time in hours; the storage 
time (90 hr being considered as maximum storage 
is included as well as the time of exposure to S—1 
sunlamp radiant energy. This graph reveals the 
stepwise pattern assumed the transmittance 
decreased on exposure to ultraviolet radiant energy, 
decreased at a different rate on storage, increased 
sharply on exposure to 2 hr of ultraviolet energy to 
the value attained preceding storage, and decreased 
on further exposure to ultraviolet energy, ete 

In order that a study of the rates of increased 
absorbance of both the pure film and the monomer- 
containing film for both the exposures to ultraviolet 
energy and the storage treatments following them 
could be made, plots of the corrected absorbance, 
A,, at various wavelengths versus hours of exposure 
to ultraviolet radiant energy were made. Figure 16 
is such a plot for 340-my wavelength. The dashed 
lines are drawn through the readings made after 13 
days of storage following the particular exposure to 
ultraviolet energy. The plots for the pure film are 
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appear to deviate from linear behavior following 
the 200-hr exposure. With up to 200-hr treatme! 
the plots for the monomer-containing film ha 
steeper slopes than those for the pure film, but 
longer exposures the slopes seem to be equal 
curves being practically parallel to each other 
The slopes of these lines plotted for numerow 
wavelengths were calculated and gave the rates 0 
increased absorbance with treatments of the film 
for the various wavelengths. These rates a 
presented in the first four columns of table 2 
an attempt to estimate the rate of increase of absor' 
ance due to storage alone, the rate due to the expe 
sure to ultraviolet energy was subtracted from th 


linear but the plots for the monomer-containing 
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8. Sources of Characteristic Absorption 
Bands in a Degraded Polystyrene 


X3M (0.0066 INCH) AFTER EXPOSUPE TO S4 


x3 (0.0066 INCH) AFTER EXPOSURE TO S$ j 
Ra ge hat 8.1. Benzaldehyde 
STORAGE ALONE OW K3A : 
It is apparent that exposure to ultraviolet radiant 


energy leads to the alteration of the polymer so that 
it absorbs in a characteristic manner, and that 
subsequent storage leads to further alteration which 
causes absorptance in a different characteristic man- 
ner. It was also noted that the absorption band 
that developed on storage was similar to the “impu- 
rity” band observed in the impure polymers. This 
similarity is observed by comparing the curve of 
impure sample X in figure 6 with the curve of pure 
polymer exposed to 155 hr S-1 plus storage shown 
in figure 12. 

Because infrared spectrophotometri studies of 
polysty rene films indicated the presence of carbonyl 
and hydroxyl groups after treatment with ultra 
violet radiant energy, and previous mass spectro 
metric measurements of gaseous products formed on 





320 degradation [3] showed some formaldehyde and 


340 
CEELSHSM, CLAREERENS benzaldehyde, benzaldehyde ® was examined in the 


of exposures on absorbance of polystyrene | ultraviolet region It gave the curves shown in 
films figure 18 for 0.1-. 1.0-. and 2 O-g liter solutions in 
chloroform. On addition of 1.6 percent of benzalde- 
ulated from the measurements made follow- | hyde to an 18-g/liter solution of pure polymer X3 
storage periods. Columns 5 and 6 of table | in’ chloroform. the spectrum obtained = (fig 19 
these rates of increased absorbance due to showed a resemblance in shape to the observed 
alone for the various wavelengths curves obtained on storage alter exposure of the 
ie rates of increased absorbance listed in | polystyrene films to ultraviolet energy However, 
ns 1, 2, 5, and 6 of table 2 are plotted versus 
negth of absorbance (fig. 17), the difference 
n the abserption curves obtained from films 
to ultraviolet energy and films subsequently 
n the dark is emphasized. The greatest rate 
e in the absorption curves obtained from 


; 


to ultraviolet energy alone is seen to be 
shorter wavelengths. The maxima in the 
ver curves indicate that 5340 mu is the wave- 


of greatest rate of change of absorbance on oF 
} | 
of the films. 
; g | 
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0005 0010 * The benzaldehyde was not purified but contained oxidation products, and 


was, therefore, considered to be more representative of the material possibly 


present in the polymer 
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the absorptance shown in figure 19 is not es broad 
and not as intense at 340 my as that which develops 
Also, there is an intense 
below 300 mu: causes little change at these 
shorter wavelengths (fig. 17 The A * of benz- 
aldehyde is 9.67, whereas its A is 0.3: therefore, 
the presence of an unreasonably high concentration 
of benzaldehyde in the sample would be required to 
give the long wavelength absorption Hence, it is 
not plausible to assume that the absorption pattern 
caused by exposure to ultraviolet can be 
attributed to the 


on storage absorptance 


storage 


ene g\ 


presence of benzaldehyde alone 


8.2. Conjugated Carbonyl-Containing Compounds 


Other compounds containing a more extended svs- 
tem of conjugation including the carbonyl group 
were therefore investigated in order that one having 
an extinction coefficient sufficient to explain either 
tvpe of absorption might be found. The APIRP 
collection of spectra [11] and the Friedel and Orchin 
collection [12] show numerous compounds which ex- 
hibit absorbance in these regions. ‘Those compounds 
exhibiting absorption bands similar in shape to those 
under discussion usually contain one or the other of 
the following structures 


he 
a 


*Kespecific extinction coefficient at indicated 
absorbance for a concentration of 1 g/liter 


avelength in millimicrons= 


Some examples are: benzoin, benzophenone, 
oxide, methyl vinyl ketone, crotonaldehyde, 
acetone, and p-toluyl acetone. However. 
tinction coefficients at the desired waveleng 
again quite small. The latter two compoun 
appreciable extinction coefficients, but the ma 
the absorption bands are shifted to 310 and 
respective ly. 

While benzaldehyde and its oxidation produ 
do not have high enough extinction coeffici: 
account for the observed in th 
stvrene films, it is likely that related struct 
as those shown below might be present in the sys' 
and their oxidation products might 
sufficient conjugation to give very 
coefficients 


absorbs ance 


easily 
large extinct 


SUBSTITUTED BENZOIN BENZIL 
There is also the possibility of a reaction similar to 
aldol condensation, once carbonyl is formed 


0 : 
" 
rr ee 


Fo hy me 
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This condensation reaction results in a final prod 
with a high degree of conjugation. Both pure al 
monomer-containing films may undergo this typ 
reaction on exposure to ultraviolet radiant energ 
once the carbonyl! is formed, 


bonyl compounds of lower molecular weight, such 4 
formaldehyde and benzaldehyde, thus 
greater mobility for reaction than might be expect 
from larger carbonyl-containing compounds. 

The yellow color evidenced by partially oxidiz 








benzaldehyde is caused by large amounts of a ma 
terial of a high degree of conjugation in the oxidatio 
Yellowing of polystyrene films can be e 


products. 
plained in a similar manner if structures like 
above are postulated. 


tioned previously. 
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Monomer styrene bias 
a reactive site for oxidation and reacts to form car- 


affording 


Infrared studies confirm tha! 
carbonyl structures do exist in polystyrene after ¢ 
posure to ultraviolet radiant energy in air as mel 
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showed no change during that time. 
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3. Stable Free Radical Compounds 


ossible that the absorbance of the film’ im- 
following exposure to ultraviolet radiant 
the absorbance of a relatively stable free 

hose structure can shift to a higher order of 

on on storage in the dark, only to be reac- 

v ultraviolet radiant energy at a later time. 

sible mechanisms might be of the types: 


On 
OARK I 
ay Rane C 
LIGHT 


SUBSTITUTED BENZOIN 
0 o 
w as 

npounds would be relatively nonvolatile at 
wratures of the exposure and storage condi- 
sed in this work. Infrared spectra of the films 
ost-radiated state were nearly identical with 
1 of the films made immediately following an 
period, Therefore, no change in carbonyl 
oxvl content between the two states could 
ected; but such a change might be too small 

easurable by this method. 
action occurring in the dark is net contin- 
sly progressive as is the reaction in ultraviolet 
ant energy in the sense that most of the storage 
tion products after any one exposure are 
d in the first 75 hr of storage and very little 
s formed even after as many as 64 days of 
ug The reaction in the dark does not seem to 
the rate of degradation by ultraviolet radiant 
but the reaction in ultraviolet governs the 
tion in the dark. In the first place, the reaction 
he dark requires a previous S-1 sunlamp exposure 
to oceur. Specimens of the pure film and of 
monomer-containing film were stored for 4 
wniths in the dark at room temperature with no 
eatment, and their ultraviolet transmittance 
Second, the 
mount of storage product formed at any one time is 
roportional to the total number of hours of exposure 





BENZIL 


to ultraviolet radiant energy received at that time. 
Even though an exposure to ultraviolet energy im- 


i 


heract tl 


nediately following a storage period seems to coun- 
ie effect of storage, the following storage 
riod shows an increase in the development of 
lorage product over the previous period due to the 


hdditional exposure to ultraviolet radiant energy. 


hese facts seem to support the assumption of the 
ormation on exposure to ultraviolet radiant energy 


s 
of relatively stable free radicals which shift to struc- 
tures of more conjugation on storage; the latter are, 
in turn, reactivated and more free radicals formed by 
further exposure to ultraviolet radiant energy. 

These “stable’’ free radicals probably can exist 
because the polymer is exposed and measured in the 
solid state. Evidence for the presence of long-lived 
free radicals in materials of rigid or crystalline 
structure is increasing day by day [14].”. The lack 
of mobility of polymer chains in a material which is 
held below its glass-transition temperature * is the 
most logical explanation for entrapment of free 
radicals on the end of polymer chains. Termination 
reactions of such radicals are extremely slow due to 
such low mobility, and the slow diffusion of oxygen 
and other gaseous materials into such a mass would 
also severely limit the rate of reaction of these 
materials with the active parts of the polymer in 
the interior of the sample. 

The presence of styrene monomer in the polymer 
apparently augments the degradation of the polymer 
chain, as well as producing monomer oxidation prod- 
ucts, such as formaldehyde and benzaldehyde 
through free radical reactions [15, 16]. In the case 
studied here, the addition of stvrene monomer to the 
polymer has resulted in degradation of the sample, 
followed by ultraviolet transmission measure- 
ments, at a rate almost double that of the pure 
polymer (see table 2 and fig. 17 Infrared measure- 
ments made on the films after the total exposure of 
290 hr to ultraviolet radiant energy from an S-1 
sunlamp showed a very appreciable difference be- 
tween the two films; the carbonyl and hydroxyl 
absorption bands of the monomer-containing film 
are almost double those of the pure film. Figure 20 
illustrates these differences well as the other 
spectral changes which occurred. The carbonyl 
group absorptance is shown in the region of 5.7 to 
6.0uand that of the hydroxyl in the region of 2.85,. 
The spectrum of each film after 290 hr of S-1 sun- 
lamp radiation is shown in this figure along with a 
spectrum of an untreated film which is common to 
both. This comparison is possible because the 
monomer that was added and remained in the one 
film was insufficient to cause any change in the 
infrared spectrum of the untreated polymer, although 
sufficient to effect an appreciable change in the 
degradation. The increase in degradation caused 
by the addition of monomer is not noticeable in the 
rates computed from the ultraviolet transmission 
measurements after approximately 200 hr of expo- 
sure to S-1 sunlamp radiant energy. This change 
in rate indicates that the added monomer had 
probably completely reacted by that time and thence 
the degradation of this film proceeded as in the pure 
film. 


as 


as 


T Such evidence was summarized by William O. Baker of the Bell Telephone 
Laboratories, Murray Hill, N. J., in a lecture at the Polytechnic Institute of 
Brooklyn on January 17, 1953, entitled ““Entrapped Polymer Radicals.” 

* Sometimes called second-order transition temperature 
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8.4. Oxidation Products 


The mechanisms that have been postulated for the 
degradation of polystyrene in ultraviolet radiant 
energy, however, are probably not all the possible 
most likely not a 
It has been shown to 


reactions because polystyren is 
linear hydrocarbon polyme 
contain oxygen in the chain [17]. Portions of cata- 
lysts when used, are also part of the chain. All these 
factors offer weak links for the initiation of degrada- 
tion and also can alter the products formed on ex- 
posure to ultraviolet radiant energy, depending on 
their type and number. The presence of a multi- 
component system is borne out in part by the fact 
that the absorption curve in the region 290 to 400 
mu for the exposed samples does not exhibit a sharp 
absorption peak, but rather a broad band, indicating 
a mixture of compounds absorbing in this general 
region 

Jjover [1S] has 
affecting the ultraviolet of a degraded 
halogen-containing copolymer. Those factors ap- 
plicable to the explanation of the absorbances noted 
in the present of degraded polystyrene in- 
clude an inductive effect of the conjugated systems 
on each other, the influence of the presence of a 
mixture of isomers, the coexistence of many degrees 
of conjugation, and the positional isomerization of 
double bonds. However, because of the presence of 
the phenyl groups and the oxidized structures in 
polystyrene, a long linear conjugated system is not 
required to produce the absorbance noted above 
300 my as is necessary in aliphatic polymers in the 
absence of air. Boyer also states that in the case 
of exposure to light, it will be those regions of the 
chain that have already developed unsaturation 
and conjugation that will absorp radiant energy 
most strongly and over a progressively widening 
spectral band. Although the oxidation and produc- 
tion of unsaturation may be a random process ini- 


described many of the factors 
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tially, the localization of absorbed energy proba 
then promotes unsaturation in adjacent monor 
units leading to the formation of further cor 
gated systems. 

The most reasonable approach to underst tain 
the steps of the reaction in ultraviolet radiant ener 
would be to study more simple parent compow 
than polystvrene. Styrene monomer, cumene 
benzaldehyde appear to be the most likely ¢o 
pounds for initial study of oxidation mechanisns 
promoted by ultraviolet radiant energy. Investig 
tions of the reactions occurring in the dark follow 
exposure to radiant energy should be studied w 
pure materials of definitely known structure wl 
are at the temperatures used. Such co 
pounds are typified py n-dotriacontane and a sim 
branched saturated hydrocarbon. Materials whos 
glass-transition temperatures are above room ten- "hy, 
perature would be desirable at the start. This 
factor would increase the possibility of obtaining 
entrapped free radicals in the mass of the mater 


solids 


9. Effect of Thermal Energy on un 
Transmittance el 

neue 

The initiation of degradation in these experiments ve 
was apparently induced by the ultraviolet radian f™....), 
energy and not by the approximately 60° C temper 7), 
ture at which the exposures were conducted. |) pe 
has been pointed out previously that 2 hr of hei)... 
at 60° C in air following a storage period es th 
which absorption occurred will not eliminate # salad 
of the absorption band formed during stor Hi. in, 
Heating at this temperature in the absence of ult... | 
violet radiant energy will not cause degradation (i pure 
the polystyrene such as that caused by expos) 
to ultraviolet radiant energy at 60° C. wanles 
Specimens of the pure film and of the monomer rified 


containing film which had been stored for 4 montis 
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irk at room temperature were heated in an 
at 60 to 63° C for a total of 290 hr. Ultra- 
insmission measurements after intermittent 
| storage periods showed no long wave- 
sorption band and almost no change at all 
yver-all spectrum. The films also did not 
ny discoloration. Previous infrared meas- 
s on films subjected to 100° C for 200 hr [2] 
how any production of carbonyl or hydroxy! 
m bands or any other spectral change. 


10. Summary and Conclusions 


wolysty rene film exposed to 290 hr of ultra- 
diant energy from an S—1 stnlamp at 60° C 
as degraded as shown by ultra- 
ransmittance. The increased 
with exposure time and the greatest rate of 
was calculated to be at wavelengths shorter 
) mu 
e of pure polystyrene film in the dark in 
oom temperature following exposure to an 
amp caused the development of a definite 
bsorption band with a maximum at 340 mug 
nt exposure to ultraviolet radiant energy 
he elimination of the absorption band de- 
on storage 
ormation of the post-radtation products was 
| by the S-1 radiant energy; storage without 
s ultraviolet exposure caused no change in 
iltraviolet spectrum. The formation of the 
products, however, did not alter the rate of 
tion by ultraviolet radiant energy. 
posure of the polymer film to heat at-60° C in 
290 hr with intermittent 
n the ultraviolet spectrum, the degradation 
by films exposed to S—1 radiant energy at this 
ture being induced by the S—1 radiant energy 
by the accompanying heat. 
ddition of monomer to pure polymer did not 
the pattern of degradation, as shown 
traviolet spectrophotometry, but did cause the 
degradation to occur at an increased rate; 
weelerated rate held until all the monomer was 
tly spent and thence the reaction proceeded 
r pure polymer. 
Reactions of the type that may account for the 
rplion curves noted on exposure to ultraviolet 
ant energy and on subsequent storage are 
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a The shape of the ultraviolet absorption curves of 
ed ‘fe purified polymer after exposure to ultraviolet 
of hes lant energy and subsequent storage closely resem- 
roe Bs the shape of the initial spectra obtained on com- 
odio real polystyrene samples, which indicates that 
od © impurities in the commercial products are of the 


ime chemical types as those produced on exposure 

polymer to ultraviolet radiant energy and | 
ubsequent storage. In most cases the commercial | 
imples absorbed more strongly than the degraded | 
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